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O ■ Abstract. The origin of prompt emission in GRBs is not yet well understood. The simplest and most popular 
> . model is Synchrotron Self-Compton (SSC) emission produced by internal shocks inside an ultra-relativistic 
\ jet. However, recent observations of a delayed high energy component by the Fermi-LAT instrument have 
' encouraged alternative models. Here we use a recently developed formulation of relativistic shocks for GRBs 
' to simulate light curves and spectra of synchrotron and self-Compton emissions in the framework of internal 
' shock model. This model takes into account the evolution of quantities such as densities of colliding shells, 
' and fraction of kinetic energy transferred to electrons and to induced magnetic field. We also extend this 
Q\ ' formulation by considering the presence of a precessing external magnetic field. These simulations are very 
. realistic and present significant improvement with respect to previous phenomenological GRB simulations. 
They reproduce light curves of separated peaks of real GRBs and variety of spectral slopes at > Epeak 
observed by the Fermi-LAT instrument. The high energy emission can be explained by synchrotron emission 
and a subdominant contribution from inverse Compton. We also suggest an explanation for extended tail 
emission and relate it to the screening of the magnetic field and/or trapping of accelerated electrons in the 
electromagnetic energy structure of the plasma in the shock front. Spectral slopes of simulated bursts at 
E <C Epeak are consistent with theoretical prediction and at £^ < Ep(,ak can be flatter if the spectrum of 
5^ \ electrons is roughly flat or has a shallow slope. The observed flat spectra at soft gamma-ray and hard x-ray 
bands is the evidence that there is a signiflcant contribution aX E < Ep(,ak from lower Lorentz factor wing 
of electron distribution which have a roughly random acceleration rather than being thermal. This means 
that the state of matter in the jet at the time of ejection is most probably nonthermal. As for the effect 
of a precessing external magnetic held, we show that due fast variation of other quantities, its signature in 
the Power Distribution Spectrum (PDS) is significantly suppressed and only when the duration of the burst 
is few times longer than the oscillation period it can be detected, otherwise either it is confused with the 
Poisson noise or with intrinsic variations of the emission. Therefore, low significant oscillations observed in 
the PDS of GRB 090709a are most probably due to a precessing magnetic field. 
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1 Introduction 

Progenitors of long GRBs are believed to be massive stars. At the end of their short lifetime of few hundred 
million years their cores collapse to neutron stars or black holes and their envelops are ejected making 
supernovae (SN) of type lb or Ic. Under special conditions that we ignore they make also gamma-ray 
bursts. However, the absence of a SN signature in some nearby GRBs such as GRB 060614 [36] is probably 
the evidence that sometimes the ejected envelop - if any - is very thin and significant amount of energy is 
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ejected as an ultra-relativistic jet which produces a GRB. We ignore details of necessary conditions for the 
occurrence of each case. 

The simplest model for the origin of GRBs is the synchrotron emission produced by collision between plasma 
shells inside an ultra relativistic jet [100, 101]. This model is a variant of fireball model which considers 
that a flash of gamma-ray is emitted when a fireball of becomes optically thin [41, 85, 91]. However, 
non-thermal spectra of gamma-ray emission of GRBs as well as broad band observations of afterglows prefer 
a baryonic fireball or a jet.^ In the latter case, dissipation of the jet's kinetic energy by collisions and shocks 
between matter shells inside the jet and their synchrotron self-Compton emission is the origin of the prompt 
emission, and the collision of the jet with the ISM produces the afterglow. In this framework a pre-burst 
tight interaction between constituent i.e. a fireball is not necessary and the main mechanism for transfer of 
energy to photons is the shock, see e.g. [97] for a review. An alternative model for the prompt emission 
is a Poynting fiow produced by energy extraction from a rotating black-hole through a Blandford-Znajek 
process [12] or from the strong magnetic field of a neutron star/magnetar [59]. In this case, the dissipation 
of energy by field lines reconnection and acceleration of electrons during this process is the origin of GRB 
prompt emission. We should also remind that although in modern GRB models the nature of the ejecta 
diverges significantly from a fireball as a dense and initially opaque concentration of energy in a small region, 
many GRB literature use the word fireball as synonym to jet/ejecta/outflow. This leads to confusion, thus 
in this work the expression fireball is used only in the context of what is called standard fireball model. 

None of the models mentioned above are flawless. In particular, since the launch of the Fermi satellite which 
is able to detect GRBs at high energies inaccessible to previous space missions e.g. BATSE and Swift, a 
number of issues and doubts have been raised against the most popular GRB models including internal 
shock-SSC, Poynting flow, and standard fireball model. In [137] authors have reviewed these problems. For 
instance, the absence of expected thermal emission from the fireball photosphere with a temperature ~ 10^ 
eV [20, 46, 86, 94, 95, 102, 119, 136] is the main drawback of both leptonic and baryonic fireball models. The 
broad band observed spectra of GRBs up to 0(100) GeV is very different from a thermal spectrum. The 
main problem of the Poynting flow model is its difficulty to explain the observed fast variations in prompt 
emission of all GRBs up to shortest time resolution of available instruments. The fiow of the electromagnetic 
energy in this model is quite smooth and probability of field line reconnection per unit time is also small. 
Moreover, simulations show that a high magnetic to matter energy at the base of the fiow becomes only 
a moderate value when the fiow arrives at a distance of ~ 10^^ cm where the formation of instabilities 
dissipates the electromagnetic energy [75, 117]. Therefore, such a Poynting fiow does not transfer the initial 
energy to a GRB as efficiently as the original suggestions expected. Internal shocks-SSC model is not without 
shortcomings either. For instance, the expected slope of the spectrum below the peak energy is mush steeper 
than what is observed [104], and doubt about this models was raised [99, 116]. Another issue for which none 
of these models has a simple solution is the observation of a tail emission in high energies E > 100 MeV. It 
fades much slower than the fiux of low energy E < 100 MeV photons[31, 32, 39, 76, 77, 96], and looks like 
to be a separate component because in some bursts e.g. GRB 090902B it has asynchronous variations with 
respect to low energy bands. 

Various explanations or solutions have been suggested for each of these problems. For instance, the absence 
of a thermal component in baryonic fireballs may be explained if the ejecta is neutron rich [65]. The decay 
of neutrons when the fireball becomes optically thin produces high energy electrons and neutrinos, and a 
significant fraction of energy is dissipated as GeV neutrinos. Moreover, Compton scattering of photons with 
electrons reduces or smears the expected thermal spectrum [92, 93]. A number of variations of Poynting 
model are also proposed [23, 63, 137] to overcome the shortages of the original suggestion. However, it 
does not seem that either of these models can reproduce the featureless spectrum of GRBs [26, 39] without 
fine-tuning. 

By contrast, it seems that internal shock-synchrotron self-Compton model can more easily explain all the 

^In leptonic fireballs a long-lasting afterglow can exist only if leptons or gamma-ray photons deposit part of their energy in 
the surrounding baryonic material. Thus a special setup of surroundings is necessary. No signature of absorption of gamma-ray 
in GRBs is observed. Therefore, observations disfavor this type of models. 



- 2 - 



issues raised here [13, 22]. For instance, in the framework of this model, the ejecta does need to be a fireball 
with a large internal energy density because, as mentioned above, the main source of energy to be transferred 
to radiation is kinematic (see [97] and references therein). Therefore, matter inside the jet can be cool [13]. 

In the context of original or tightly coupled baryonic fireball models, the role of the photosphere is 
crucial because GRB gamma-ray emission is considered to be due to the simultaneous release of trapped 
photons when the fireball becomes optically thin. Even if the ejecta/jet is a fireball, the absence of a thermal 
emission can be explained if internal shocks occur inside the photosphere. In this case, chaotic collision, 
acceleration of particles, expansion, and finally emission of synchrotron photons can completely smear the 
initial thermal state of the ejecta. The radius of the photosphere of a standard fireball is estimated to be 
i?p/j ~ 6 X 10i^(L/10^2 gj.g sec-^)C"^(r/100)~2 cm [20, 26, 86, 102], where L is the total luminosity of the 
baryonic outflow of the fireball [26], C, = Lin/(?M is its initial dimensionless entropy, and F is the Lorentz 
factor of fireball [20]. From the Fermi [7] and Swift [34] observations a Lorentz factor of F > 1000. If total 
outfiow luminosity is L ]^q53-55 gj,g ^^^^ ^ \ _ j^qq^ photospheric radius Rph is in range 10^*^"^^ cm [20]. If 
thermal radiation is the only or the dominant contributor to observed radiation, it must have a temperature 
~ 1 MeV. Under some conditions (see next paragraph) internal shocks can occur inside this radius. Note 
that in internal shock+SSC models the ejecta is usually considered to be cold and there is no issue with the 
absence of a thermal component. In these models shocks usually occur outside the photospheric radius and 
the photospheric emission is weak and at low energies. 

The distance of internal shocks from the central engine is estimated by the time necessary for a fast shell 
with Lorentz factor Fj to catch and collide with a slower shell with Lorentz factor F^ ejected At sec earlier: 

2cAtF^F^ 

RcoU = p2 _ V2 - ^^^^r' r/»F, (1.1) 

All quantities in (1.1) are in the rest frame of the central source. GRB literature usually assume At ~ 0.1 sec 
as the typical time separation between shells - density anisotropies in the fiow. This leads to Rcoii ~ 10^^ cm 
i.e. the same order as the photosphere radius of a standard fireball for some choice of parameters - specially 
when a low F is assumed. However, the choice of At ~ 0.1 sec is arbitrary because in both long and short 
bursts we see variations at much smaller time scales. In fact, the whole duration of some short bursts is less 
or comparable to this value. On the other hand, simulations show that collimated jets come from the most 
inner part of the accretion disk. The size of neutron stars and stellar-size black holes is few kilometers. Thus, 
the innermost radius of the accretion disk must be of the same order. Assuming that spatial separation 
between density anisotropies of the disk is comparable with this inner radii, the natural variation time scale 
for relativistic accretion disk is ~ {Rdisk/'^ km) x 10~^ sec. This leads to Rcoii ~ lO^''"^^ cm for progenitors 
of mass ~ 2 — 200 solar mass. Even if we consider the orbital period of spinning material around black 
holes with masses a few times Mq, the Keplerian approximate period for the innermost quasi-stable orbit 
r = (47rGM/c^)^/^ is of order ~ 10~^ sec. This leads to RcoU ~ 10^^ consistent with the upper limit of 
the range given above. Besides, simulations of magneto-rotational instabilities in accretion disks show the 
presence of a "channel" type radial flow both inward and outward [58]. This radial motion is necessary for 
reducing the viscosity of the accreted material which suppresses the conductivity of the fiow and thereby 
the formation of a strong magnetic field and a jet. Therefore, Rcoii range obtained by assuming the radius 
as the distance scale between anisotropies is reliable. Detail modeling of some bursts with good data (see 
e.g. [83, 140]) also confirms these estimations^. 

If the ejecta is a standard fireball, short distances suggested here for internal shocks mean that their initiation 
inside the photosphere is not ruled out if the input entropy C, and the Lorentz factor are small or if At is 
very small. Nonetheless, following the internal shocks, acceleration of electrons, formation of a random 
magnetic field, synchrotron/self-Compton scattering, and sudden blow up of the ejecta, the initial strong 
coupling and thermal distribution of particles will be interrupted and their memory in the distribution of 

■^We should remind that in GRBs variabihty is observed in all scales. This means that the accretion disk and thereby the 
jet don't have a special variation distance scale and are roughly scaleless. This ensure that the process of internal shocks is not 
limited to small scales and GRB emission as a whole is dominantly due to this process. 
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high energy electrons will be in large extend lost [111], see sections 4.4 and 4.4.1 for more details. Therefore 
in general during the whole or large fraction of burst time no photosphere is present and photons are free to 
propagate. We should also mention that according to PIC simulations despite metamorphoses mentioned 
above, a dominant low energy thermal population of electrons seems to persist, at least in plasma. On 
the other hand, internal shock models do not need a fireball, thus issues related to shocks inside or outside 
the photosphere do not arise. This scenario is consistent with non-thermal spectrum of prompt emission 
in high energy bands. Other issues of internal shock-SSC model such as small efficiency, steep slope of the 
spectrum at low energies which is not observed, etc. are related to the micro-physics of the shock and can be 
realistically addressed only when more realistic PIC simulations become possible. If internal shocks occur 
at short distances from the central source, the densities of the shells and surrounding are expected to be 
large and optical depth small. We address this issue in Sec. 4.2 

In order to better understand internal shock-SSC model in a previous work [143] we have developed a simple 
and consistent formulation of ultra-relativistic shocks. This formalism and corresponding simulations are 
essentially similar to work by other authors specially [13, 22]. The latter authors have emphasized on the 
phenomenological aspects of the formulation and given a qualitative view on the role of various quantities 
on observable properties of GRBs. For this reason these topics have not been discussed in [143] or here. On 
the other hand, we have tried to discuss in details mathematical aspects of the problem as much as our aim 
to keep the formulation simple allowed. Moreover, this formulation includes time variations of micro-physics 
parameters such as densities and fraction of energy transferred to accelerated electrons and to the induced 
magnetic field in the shock which in the context of fast varying prompt emission are crucial. The advantage 
of our method is that it can be relatively easily extended and time variation of other quantities can be 
included. In addition, giving the fact that at least 3 moving reference frames are involved in the modeling of 
relativistic shocks, a detailed formulation reduces the risk of confusion for the reader. For each simulation 
we present light curves and spectra such that they can be directly compared with observations. Due to much 
larger parameter space in these simulations with respect to the modeling of [13], we leave its exploration 
to a future work. Nonetheless, we have included results for bursts with different durations, brightness, 
and hardness to show that practically all the features seen in real GRBs can be reproduced. Notably, the 
variety of spectral slopes which has been recently claimed as to be unexplainable with the simple internal 
shocks-SSC model are obtained. 

Preliminary simulations of both prompt and afterglow according to this formulation [144] have shown that 
they well reproduce the behaviour of GRBs in various intervals of their evolution. Here, we improve those 
results by simulating the whole prompt emission of GRBs. In addition, we extend the formalism of [143] 
to include an external magnetic field, and we study its effect on the prompt emission of simulated GRBs. 
Here by external field we mean not induced by shocks. It can be of one or a combination of the following 
origins: 

1. A Poynting flow i.e. a slowly varying electromagnetic wave; 

2. A magnetic field frozen in the jet's plasma i.e. the magnetic field in a magnetohydrodynamical flow; 

3. Magnetic field of the central source in case of a proto- neutron star, magnetar or accretion disk around 
a black hole; 

4. Remnant of the dynamo field in the star's envelop. 

Considering the distance to central source obtained earlier, shocks can be surrounded by this envelop if 
it is not yet ejected. In the case 1 precession of the field seems unlikely. In case 2 the field precesses 
only if the plasma has a helical twister- like movement. In cases 3 and 4 the magnetic field can have a 
precession. In case 3 the precession period is expected to be of millisecond scale. If the outer layers of 
the star continue their differential rotation before being completely disrupted and ejected, their dynamo 
magnetic field (case 4) should have a rotation period of few seconds. We do not specify the type of the 
field in the simulations and consider the general case of a precessing field. A non-precessing field can be 
assumed to have a precession period much longer than prompt emission of simulated GRBs. The simple 



phenomenological model used for simulations here is not sensitive to these details. Thus we simply assume 
that the magnetic field influencing accelerated electrons has 2 components: one is induced by the shock and 
is parametrized as usual by fraction of energy of falling ions transferred to electrons, see equation (3.17) in 
the next section, and the other component is assumed to be independent of the shock and a function of the 
distance from central source only. 

We have also calculated the distribution of Compton scattered photons by electrons for these simulations. 
We show that its contribution is very small. This is consistent with observations and is another evidence 
of the validity of internal shock-synchrotron model. Compton scattering between protons and photons, 
photon-photon annihilation, and Compton self-absorption are not considered. It is known that Compton 
emission extends longer in time than synchrotron, but we show that it cannot explain observations by the 
Fermi-LAT [32, 96]. On the other hand, based on the results of the state-of-art Particle-In-Cell (PIC) 
simulations we related the high energy tail to the structure of induced electromagnetic field in the shocks. 
It can play an important role in a delayed dissipation of most energetic electrons. The screening of the 
external magnetic field by the shock can be also involved. This probably explains why the delayed tails 
have been observed only in bright hard bursts in which, as we will show here, an external field can have 
significant contribution. However, the phenomenological model of ultra-relativistic shocks used here cannot 
address these complexities. 

We begin by briefly discussing evidence for the presence of precessing magnetic fields in astronomical sources 
and the case of GRB 090709A in Sec. 2 as an observational evidence in favour of an external magnetic field 
in GRBs. The formulation of dynamics in presence of an external field is discussed in Sec. 3. In Sec. 4 we 
present a number of simulated bursts, discuss their properties and compare them with properties of observed 
bursts. This section also includes results of the simulation of Inverse Compton contribution in light curve 
and spectrum of GRBS, as well as an explanation for the high energy tail emission. Finally Sec. 5 includes 
outlines of this work. In Appendix A an analytical expression is given for the contribution of the external 
magnetic field in the evolution of bursts kinematic. In Appendix B the projection of a precessing magnetic 
field on the shock front surface is obtained. Finally, Appendix C gives a qualitative description of how a 
magnetic field influences the formation of instabilities in charged plasmas. 

2 Precessing magnetic fields in astronomical sources 

There are a number of motivations for presence of other sources of magnetic field than shocks in the 
environment of relativistic jets. For instance, due to turbulent rotation of the accretion disk that surrounds 
the central source the disk is highly magnetized. Along with dynamical and gravitational forces, the pressure 
created by this field plays an important role in the formation and collimation of ultra-relativistic jets (see e.g. 
[53] for a review). In fact observations of astronomical sources in which long lasting jets are formed such as 
Active Galactic Nuclei (AGNs) [15] and accreting X-ray binary stars [60, 80], tell us that both mechanical 
and magnetic pressures must be involved (see e.g. [4] for a review). In all these cases jets are associated 
with a rotating and precessing accretion disk [78]. The mechanical pressure at the inner side of the disk 
compresses and pushes part of the falling material up and ejects it out to free the remnant momentum and 
angular momentum of in-falling matter. Therefore, it seems that the presence of an external magnetic field 
in GRBs environment is inevitable. 

In addition to providing acceleration and collimation, a strong magnetic field explains the puzzling and 
apparently inconsistent fluence of high energy component of the bursts observed by the Fermi-LAT instru- 
ment. Simulations presented in the next sections show that when an additional magnetic field is included, 
bursts are few times brighter and harder. This soften the efficiency issue raised by Particle In Cell (PIC) 
simulations which find that only ~ 10% of kinetic energy of the jet is transferred to electrons [111]. Although 
a larger magnetic field does not solve the problem of insufficient energy transfer from bulk of material to 
electrons, it makes synchrotron emission more efficient. The low efficiency of energy transfer means that the 
total energy released by the engine must be even larger. However, simulations by [111] are performed for a 
plasma. Due to the large mass difference between nucleons and electrons, the behaviour of a baryonic 
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plasma can be very different from a e one. Indeed PIC simulations with larger mass difference between 
particles confirm this point. In particular they show that the presence of a parallel magnetic field facilitates 
the acceleration of electrons [72]. 

Due to the fast precession of accretion disk and probably the central source itself, it is expected that during 
a coUapsar phase the external field at the position of internal shocks precesses. This can be the case even 
if the jet is aligned with the magnetic field at the time of its ejection. Moreover, if the argument given in 
the Introduction is correct and GRBs occur in a distance of ~ 10^''"^^ cm, at the time of their formation 
internal shocks are yet inside the star's external envelop, thus the flux of the field can be significant. 

2.1 A special case: GRB 090709 A 

The absence of a clear evidence of a precessing field in prompt emission of GRBs was a puzzle, but the 
situation has been changed with the observation of GRB090709a by Swift-BAT, [62], Konus-Wind [40], 
INTEGRAL-SPI-ACS [42], and Suzaku-WAM [82] instruments. Initial doubts on the identification of this 
transient as a GRB has been removed by further observations [68] including the detection of a candidate host 
galaxy for this burst [90]. The signature of an oscillating component with a period of ~ 8 sec superposed on 
two FRED-like spikes is detected by all the instruments that observed the prompt emission in gamma-ray. 
Although from further analysis of the data [16] it has been claimed that the detected oscillations in the PDS 
have only a significance of < 3a, the fact that they are detected by multiple independent instruments and 
in different energy bands is the proof that they are genuine. Oscillations are particularly significant in the 
power spectrum of light curves rather than in the differential spectrum obtained after removing a smoothed 
curve. The latter is necessarily depends on the way the smoothed component is determined and can include 
artifacts to power spectrum. 

Fig. 1 shows the power spectrum of GRB 090709A in 4 non-overlapping Swift-BAT bands. In place of mod- 
eling the noise, we have calculated the random variation of the average spectrum by adding a Gaussian noise 
with a standard deviation equal to the estimated error for each data point. We have repeated this procedure 
500 times. In Fig. 1 the plateau formed by the simulated noise shows precisely the level of white noise in 
the data. Features covered by points are not distinguishable from noise. The peak around / = 0.12 Hz is 
about one order of magnitude higher than the noise in all 4 non-overlapping energy channels. Considering 
the number of iterations, we obtain a significance of 3.5a > 99%) per channel for this feature. Therefore, 
assuming that the measurement errors in different channels are independent, the probability of coincidence 
in all 4 channels is < 10~^. We have performed the same analysis for a number of other GRBs detected 
by the Swift-BAT instrument and didn't find any other burst showing an oscillatory component in its PDS 
with such a high significance. The result of this study will be reported elsewhere. No significant evidence of 
oscillation was found in the Swift-XRT Window-Timing (WT) data beyond 150 sec after trigger [68]. This is 
consistent with our simulations which show that a significant oscillation exists only in the very early prompt 
emission because when the jet leaves the envelop of the star, the magnetic field decreases very rapidly. 

To interpret these oscillations we briefiy review their plausible sources. An oscillatory modulation of emission 
from an astronomical source can arise from a number of processes: 

• Precession of a hot spot: In the framework of internal shocks and synchrotron/inverse-Compton 

emission model such explanation does not apply to GRBs. After the ejection of the relativistic jet, 
gravitational and kinematic infiuence of the central engine on the ejecta are negligible. The exception 
could be the emission very close to a rotating black hole. In this case the relativistic ejecta could 
be dragged by the hole for sometimes after its ejection. For instance, the ejecta can have a helical 
orbit [118] before the gravitational infiuence of the hole or precession of a jet around the rotation 
axis of a Kerr black hole is stopped. In this case the effect should be energy dependent because less 
energetic photons should be less affected by the black hole drag. 

• Regular and oscillation-like density anisotropics in the ejecta Density perturbations in the ejecta 

are expected and are one of the main ingredients of the internal shock model. However, it does not 
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Figure 1. Orange (light grey): Power spectrum of GRB 090709 A in 4 non-overlapping Swift-BAT bands, from 
top left to bottom right: 15-25 keV, 25-50 kcV, 50-100 keV, 100-350 keV. We have used BAT light curves with 64 
msec resolution and binned them (bin-size = 10). Turquoise (mid-gray) points present 500 spectra of modified light 
curves by adding to each data point a random value with Gaussian distribution and a standard deviation equal to the 
estimated la error of the Swift-BAT data. Purple (dark grey) line presents a fit to obtain an average slope index of 
a ~ —1.3 of the spectrum (for an ideal bilateral-triangle shape peak a — —2). The peak around / = 0.12 Hz is clearly 
distinguishable from the noise and is consistent with the analysis of other groups. There are also some other peaks at 
higher frequencies but with much less significance. Inspection of plots with linear time scale shows that most probably 
they are not noise but the harmonics of / = 0.12 Hz oscillation. The two broad peaks at low frequency limits are due 
to the main peaks in the light curves and their first harmonic. 



seem plausible that anisotropies be able to keep their coherence for a long time to induce a regular 
modulation in the emission. Expansion of the ejecta, turbulence, and minor internal shocks should 
very quickly destroy the coherence of different bunches leaving only random variations - substructures 
- in the emission. Nonetheless, if the source of oscillations is external, for instance if they are produce 
under the influence of an oscillating magnetic field, distribution of the particles in the plasma can be 
modulated (see also Appendix app:c). 

• External magnetic field: The presence of a pre-collapse large scale strong magnetic field, in the range 
of few hundred to few kG and as large as few tens of kG in 0-stars and Ap /Bp stars respectively, has 
been confirmed by many observations (see [45, 127] for a review). Moreover, the iron core of collapsars 
progenitors are expected to have a fast rotation with a period of ~ 1 — 8 sec [120]. The combination 
of an external magnetic field and rotation not only helps but perhaps is indispensable [9, 69] for the 
formation of an ultra-relativistic jet [14, 52] and a GRB. 

Therefore, it seems that a rotating magnetic field is the most plausible explanation for the cyclic modulation 
of peaks in GRB 090709a. In fact if our understanding of the structure of massive stars and supernova 
explosion is correct, the effect of magnetic field of the central engine should have been seen in most GRBs. 



In Sec. 4 we show that due to fast variation of quantities that control the luminosity of a burst, coherent 
oscillation of quantities such as magnetic field is smeared or destroyed and even when they are present in 
the light curves, only under special conditions they can be distinguished from from intrinsic fluctuations of 
the emission. 

3 Synchrotron emission by relativistic shocks in presence of an external magnetic field 

In this section after a quick review of the simplified relativistic shock and synchrotron model discussed 
in [143], we reformulate the model in presence of an external magnetic field^. 

A shock is formed when two shells of material with different densities collide. The discontinuity in physical 
properties of the shells exist only at the beginning of the collision, thus is considered as an initial condition 
(see e.g. [2]). During collision one can distinguish two shocked zone in the opposite side of the initial 
discontinuity. For an observer at rest with respect to the shock front - which is no longer a discontinuity but 
the place of a sharp variation of density and other physical properties - some of the particles are accelerated 
and ejected to downstream. On the upstream side by contrast particles are decelerated. This process 
induces turbulence and creates coherent electric and transversal magnetic fields which accelerate electrons 
by Fermi processes [11, 72, 111]. These electrons move back and forth across the shocked zone. Their 
energy is dissipated mainly by synchrotron emission in places where the induced magnetic field is strong 
and transversal. If the velocity of massive particles - presumably baryons - are reduced to relativistic sound 
speed in the upstream, a secondary reverse shock front which moves in the opposite direction - with respect 
to the main forward shock would be formed [133]. 

Although it was expected that the difference between the dominant synchrotron frequency and time evolu- 
tion of light curves from forward and reverse shocks make their separation possible, multi-band and early 
observations of GRBs showed the contrary. In fact since the first GRBs detected by Swift there was a 
significant effort by observers to detect what is called an optical flash with steep slope characteristic of 
the emission from reverse shock according to some theoretical models [74, 107]. Despite interpretation of 
some bright early optical afterglows as being the prompt reverse shock, none of them have the necessary 
properties - the rising slope - to be indisputably interpreted as a reverse shock peak. In all such cases either 
the observed peak in optical coincides with flares in X-ray and/or gamma-ray, or can be interpreted as the 
passage of the synchrotron characteristic frequency of the continuum to optical bands followed by a fast 
decline [79]. A good example is GRB 061121 [87] in which the burst was detected during its precursor and 
optical data are available during its main peak. Simulations of shocks (including ours in this work) reproduce 
this continuum-like peak in low-energy bands few tens of seconds after the prompt peak in gamma-ray [67]. 
It is produced by forward shocks. Simulations in [67] also show the reverse shock peak occurs before forward 
shock peak. However, such a peak has not been observed in any bursts with optical observations before a 
major peak such as GRB 050820a [89, 126], GRB 061121 [87], and GRB 070721B [142]. Therefore, based on 
these observations one can conclude that in the shocks which make the prompt emission of GRBs despite 
energy dissipation in the shock, the velocity of particles in the jet is not reduced down to relativistic sound 
speed such that a reverse shock be produced [133] or the reverse shock is too weak to be detectable, or the 
delay between arrival time of forward and reverseshocks emission to observer is too short to make them 
distinguishable. As for PIC simulations which are more precise, shocks in shows very little reverse shock, 

^As we have mentioned in the introduction what we call an 'external magnetic field' can have various origins, and is not 
necessarily a Poynting flow. Moreover, relatively small fields we consider do not make the jet a 'magnetized outflow' i.e. Poynting 
energy dominant (see Fig. flg:kinetic), and a dust approximation is adequate in what concerns the kinematic of the flow. We 
also remind that this work is not meant to be a magnetohydrodynamics modeling of GRBs which is definitely a much more 
complex problem. Apart from PIC simulations some authors have used simplified magnetohydrodynamics algorithms to study 
dynamics of relativistic shocks between magnetized plasma shells without considering any radiation emission [66] . In their model 
synchrotron emission is added analytically [67]. In fact, even in absence of an external magnetic field, a correct formulation of 
the jet needs a magnetohydrodynamics modeling because most probably the jet is ionized, includes current sheets, and thereby 
a magnetic field - although it can be small and has a small coherent length. Our analytical model and simulations presented 
here are evidently much simpler, nonetheless, they allow to investigate microphysics parameters and can be used to model real 
bursts. 
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see Fig l-e,h,i in [111]. No realistic simulation of baryonic relativistic shocks is yet available, see e.g. [72], but 
one can argue that if light particles such as electrons with large e/m can be barely reflected by a relativistic 
shock, it would be even more difficult for massive particles like protons. 

In this case emissions from the prompt shock is better modeled by assuming a single emitting region that 
we call the active region. It is dominantly in this region that transfer of kinetic energy from fast shell to 
slow shell and radiation occurs. The active region is expected to be smaller than shocked region - the zone 
in which two shells are mixed. Notably, in contrast to the shocked region which increases continuously 
until two shells are coalesced or pass through each other, the width of active region first increases, arrives 
to a peak value, then decreases when the discontinuity between colliding shells changes to a continuous 
transfer region. To simplify the model further we also assume that the thickness of this emitting region is 
small, thus the propagation time of photons in this region is smaller than time resolution of this model. 
In fact for objects moving with ultra relativistic speeds with respect to an observer, time and distance are 
approximately proportional: r'{t') = I3'{t)ct' w ct'.^ Under these approximations evolving quantities only 
depend on the average distance of the active region from central engine (or a far observer). Mathematically, 
this approximation is equivalent to assuming a wave-like behaviour for dynamical quantities i.e they depend 
on r' — c/3't' rather than r' and t' separately. When /?' = const, i.e. when there is no collision or dissipation, 
this is an exact solution. In this case the solution at every point can be obtained from the solution of one 
point. 

Conservation of energy and momentum determines the evolution of this system. The velocity /?' of the 
fast shell/active region decreases due to absorption of particles from slow shell and dissipation of kinetic 
energy as radiation due to synchrotron and self-Compton interactions. After a variable change the dynamic 
equations - energy-momentum conservation equations - for the active region can be written: 

dir'-'n'Ar'i) _,f,,d{n'Ar') , , ,2. , a .'a _ , , ^ ,2 dE', 



^ + 2r'(n'Ar')) + r'\n' Ar')^ = n',ir)r- - (3.1) 



dr' ^ \ dr 

, , = . , + 2r'(n'Ar')) + r'\n' Ar')^^ = — (3.2) 

dr dr dr ATrmc^dr 

where r' is the average distance of the active region from central engine, n' is the baryon number density 
of the fast shell measured in the slow shell frame, tiq is the baryon number density of the slow shell in its 
rest frame and in general it depends on r'. Here we assume that ng(r') = NQ{r' /rQ)~'^. For the ISM or 
thin shells where density difference across them is negligible k = 0, i.e. no radial dependence. For a wind 
surrounding the central engine usually k = 2 [17] is assumed. For a thin shell/jet expanding adiabatically 
also K = 2 if we neglect the transverse expansion in the case of a jet (collimated ejecta). If the lifetime of 
the collision is short we can neglect the density change due to expansion during the collision and assume 
K = 0. Ar' is the thickness of the active region, 7' is the Lorentz factor of the fast shell with respect to the 
slow shell, /3' = 7'^ — I/7', m = rup + rUe ~ nip, E'^y is the total emitted energy, and c is the speed of 
light. The evolution of the average radius of the shell is: 

r'{t') - r'(t[,) = c / I3'{t")dt" (3.3) 
where the initial time t'^ is considered to be the beginning of the collision. 

If the thickness of the active region is not negligible, as a first approximation, we can assume that it can 
be divided to layers with different characteristics which are in a quasi-steady state with respect to each 
others. In this case equations (3.1) and (3.2) can be written separately for each layer and the emission can 
be integrated. Note also that the thickness of the shell in these equations always appears as = n'Ar', 



^Through this work quantities with a prime are measured with respect to the rest frame of the slow sheU and without prime 
with respect to a far observer at the redshift of the central engine. Parameters used for parametrization do not have a prime 
even when parametrization is in the slow shell frame. 
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i.e. the column density of the sheh. Thus in the hmit of shells with infinitesimal thicknesses we can simply 
replace n'Ar' with n'^ which is in addition a Lorentz invariant quantity.^ 

After integration of equation (3.1) we obtain the evolution equation of the column density of the active 
region: 

, Kr'li^v^,?-'' - 1) + (3 - K)/g<(/o)7^(l - P',) 

(3 - K)r'2y(l - /?') ^"^^^ 
where 7q = "{'{r'^) and /3q is the corresponding /3'. By solving this equation along with (3.2) we obtain the 
evolution equations for 7' and n(, with respect to radius or equivalently time. 



P' = -^ = c/3'-f = ^r'^^r'arn'^^ / <{lehld^e (3.5) 



The power of synchrotron radiation emitted by the active region is: 

dE' a, dKy 167r ,2. / /2^i^ 

where n'^ is the number density of accelerated charged leptons - electrons and possibly positrons - with a 
Lorentz factor 7e 7' with respect to slow shell frame. We assume that the thermal motion of charged 
particles is negligible with respect to their boost, thus, B'^ is defined as the total magnetic field component 
perpendicular to the direction of accelerated electrons. B'j_^ /Sir is the magnetic energy density and ax is 
the Thompson cross-section. 

The normalization of the electrons distribution is defined as: 

"-6(76)^76 = < (3.6) 

7m 

00 ^,'2™, „/ 



-feneilejd-ye = (3.7) 



jm ^ 

where n'^ is the number density of accelerated charged leptons and eg is the fraction of the kinetic energy 
of the falling baryons transferred to the accelerated leptons in the slow shell frame. Here we assume that 
it is equal to the number density of the slow shell accelerated to fast shell i.e. n'^ = j'tiq. Note also that 
(3.7) assumes an abrupt break in the electrons Lorentz factor distribution. This simplifying assumption has 
crucial consequences for the spectrum of synchrotron emission that we will describe in Sec. 4.4.1. 

From constraints (3.6) and (3.7) the minimum Lorentz factor of electrons jrn and the normalization of ng(7e) 
can be obtained. For a power-law distribution: 

n'M = iVe — for 7g ^ 7^ (3.8) 

V7m/ 
/ 2 /2 

Ne = — = 7 ^j^^ r 3.9 

{P - l)eej'^mpn'Q 

7m = ] (3.10) 

In the simulations presented here we also consider the case of a distribution with an exponential cutoff at 
high energies: 

\-(p+i) 

ng(7e) = Ne( — ] exp( —) , 7e ^ 7^ (3.11) 



7m / 7cttt 
te7 ^p^o 7c„i V y^^^. ) 



X ^ (3.12) 

^e = 7 ' r(a,x)=/ dyy'^-^e-y (3.13) 

7n 



^During a collision the active region is essentially part of the overlap of two shells. In a radiative collision or even when the 
exchange of energy between shells is purely elastic, this region has necessarily a smaller Lorentz factor than yet uncoUided part 
of the jet. 
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In simulations present in next sections we assume that ^m/lcut factor is constant during each regime. 
Therefore in equation (3.12) this factor is a parameter rather than a variable. From properties of incomplete 
Gamma function we find that when ^mllcut — ^ 0, the solutions for 7^ and Nf. are approximately equations 
(3.12) and (3.9). We also use a broken power-low distribution for electrons in the simulations: 



\ -(pi+i) 

^e[T^\ 7m ^le^ Icut 

\ -(P2+l) / \P1-P2 

^e(^ le>lcut, K = nA^^ 



(3.14) 



Ne = 7 ^ (3.15) 



' PimeU' 1 - (1 - anl)!- 7™_)Pi-l ^ ' 

In absence of an external magnetic field the total field B' is equal to the magnetic field induced by turbulence 
and Weibel instabilities [130] in the shocked material (the active region) Bj^^. Its direction is almost per- 
pendicular to the direction of motion of charged particles i.e. to the surface created by the shock front [72].^ 
This means that Bjj^_|_ Bj^^. The amplitude of Bj^^ is expected to be proportional to the energy density of 
particles falling into the shock front and can be parametrized as: 

^ = eBC^^'^m^n'^ (3.17) 
where is the fraction of protons energy transferred to the induced magnetic field. 

In presence of an external magnetic field the total field is B' = B-j^ -|- B^^. In this case the amplitude of the 
magnetic field in the shock front surface is: 

B'^ = B'l^ + B'l^ + 2B',^.B'^^ (3.18) 

If the external magnetic field precesses around a fixed axis, both its amplitude and its angle with the shock 
front surface S are time dependent and periodic. On the other hand, B-^^ is a circularly polarized random 
field across the shocked zone (see e.g. both ID and 3D simulations in [72]), and B'in presents the mean 
value of its amplitude. Its projection on S is also random with a uniform distribution if the ejecta has a 
cylindrical symmetry - the case we consider here. Therefore, the angle between Bj^^ and B^^ is random. 
If its random variation takes place much faster than the precession of the external field, the average scalar 
product of (Bj^ • Bgx) = 0. Therefore, for determining the total power of synchrotron emission and the 
spectrum we only need to determine the projection of the precessing field on the shock front surface. The 
details of the calculation can be found in Appendix B. 

From dynamic equations (3.1) and (3.2), the solution obtained for the column density n'^ (3.4), and the 
relation between synchrotron power and the amplitude of the magnetic field (3.5), we obtain the following 
equation for the evolution of /?': 



d 



d 



iff- - 1 + (^-)";;ro)^-'(-n) (i _ ))^. 



(3-k)(1-/30 



Ai^^r' (r'\ 2-'' Aii^^r' ( r' ^ ^'^^ 



/3'Ar'(r^) Vro/ /3'Ar'(r^) 



(3.19) 



. _ 4amgarn^Ar^(r^)e^(r^)eij(r^) 

3m2 ^"^-^"^ 

_ «"^P^TA/(r^)e^(r^)g^e^_L(r^) 



^In reality electrons have a helical movement in the direction parallel to the jet axis [24], but we neglect these details in the 
simple model presented here. 



where 77 = la^ + a_B + 2k and r\\ = 2ae + k + ax, where Oe, as, and ax are respectively power-law indices 
of Ee, eB, and external magnetic B^x defined as x oc (r/ro)"" for each of them. 

Note that the two terms in the right hand side of (3.19) have very similar forms. Although this differential 
equation is first order and can be solved easily, it is also highly nonlinear. Indeed, the solution includes a 
polynomial of order 6 of /?' and a term proportional to (1 — /3')^/^. Such an algebraic equation must be 
solved numerically at each point and is not suitable for numerical simulations. Another approach - which 
is explained in details in [143] - is a perturbative/iterative approximation. It consists of the integration of 
right hand side of (3.19) by replacing f3' with its order n — 1 approximate solution to obtain its order n. 
The solution of order zero corresponds to absence of synchrotron emission, i.e. the right hand side of (3.19) 
is set to zero. The general solution at order n has the following form: 

p _ •AA4(„^i)(rO+^iA/-(„_i)(rO 

«' = ^'•'(^0) CO noN 

f^in) 1 ^^r'^3->.-, 1^ I V ■4A4(„„i)(rO+ytiA-(„_i)(rO ^"^'^^^ 

3-ftVW'o'' ^ Ar'(r^,) 

^..)(/)^^/^^.-.. (3.3) 

^.-,(0^^/^^.-M. (3.24) 
^ ^ n'ir',)Ar'ir',)/3',% 

The advantage of this approach is that one can use the approximate analytical solution to study the effect of 
various parameters and quantities on the evolution of dynamics of the ejecta and its synchrotron emission. 
To proceed, we have to consider a model for the evolution of Ar'(r') which cannot be obtained from the 
simple formulation of relativistic shocks presented here. We consider the following phenomenological models: 



Ar' = Ar'^ 
Ar' = Ar^ 



r 



1 — exp(— J' — 



e(r' - r'o) Steady state model (3.26) 
0(r — Tq) Exponential model (3.27) 



Ar' = Aro^^^ Q{r' - r'^) dynamical model (3.28) 

Ar' = Ar^^^^ Q{r' - r'^) Power-law model (3.29) 

Ar' = ArQ exp | —S'— ) 0(r' — Tq) Exponential decay model (3.30) 

The initial width Ar'(rQ) in the first two models is zero, therefore they are suitable for description of 
initial formation of an active region in internal or external shocks. The last three models are suitable for 
describing more moderate growth or decline of the active region. The integration of (3.23) and (3.24) for 
models containing power-law terms can be obtained from the solution for dynamical model. The results 
for A4Q{r') is presented in [143] and for (3.24) can be found in Appendix A. For models containing an 
exponential term analytical calculation is too complex and only numerical integration is used in our studies. 

For the simulations described in the next section we have tried both approximated analytical solutions for 
A4Q{r') and Mo{r') - in the cases for which an analytical solution is available, and numerical integration of 
(3.23) and (3.24). We have also used the exact solution of differential equation (3.19) and its linearization 
to find an approximate solution for /?'. The results of these 3 methods differs by a few percents. Therefore, 
considering other simplifications and uncertainties in this model, it is not important which one of them is 
used. However, in practice it is easier to use the approximate solution (3.22) with numerical integration of 
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integrals because this procedure can be applied to all models for Ar'(r') and new models can be added to 
simulation code without additional analytical calculations. 



The next step is the determination of the synchrotron emission spectrum and flux. In [143] (and references 
therein) it is shown that energy (intensity) angular spectrum of synchrotron depends on i?^ the transverse 
component of the field with respect to the direction of motion of electrons only through the dependence on 
the synchrotron characteristic frequency: 

u' = — 3.31 

2cme ^ ' 

Therefore, when an external magnetic field is present, the formulation developed in [143] can be used 
without any modification. We only need to consider the total field in the determination of w^. For the sake 
of completeness we repeat here the final expression obtained for the synchrotron emission in [143]: 



dP AVSe^ , Ar 

= r 

ooduj 3tt r(r) 



/ d7e<(7e)7e~'^2/3(-)+-^(^,0 (3.32) 

" y?T). c 



where -F(w, r) includes subdominant terms and terms depending on the curvature of emission surface. 
In [143] it is argued that these terms are much smaller than the dominant term, thus we neglect them in 
the simulations. The slow shell Lorentz factor r(r) can be expressed as: 

r(r) = r/7'(r)(l + /3//3'(r)) (3.33) 

where Tj is the Lorentz factor of coalesced shells with respect to far observer at the same redshift. Note 
that synchrotron emission spectrum in (3.32) is integrated with respect to to angular distribution. The 
description of angular dependent spectrum can be found in [143], equations (56) and (60) for slow shell and 
observer rest frame respectively. We use (3.32) in our numerical simulation rather than angular dependent 
spectrum because it makes the comparison with observations easier. However one should remember that 
the total emission energy obtained from this formula is Eiso and one has to take into account the effect of 
relativistic beaming of the emission in the rest frame of the observer if one wants to compare the energy 
output with the initial kinetic energy of the ejecta, see Sec. ?? 



4 Simulation of light curves and spectra of GRBs in presence of an external magnetic 
field 

In this section we present light curves and spectrum of simulated GRBs according to the simplified relativistic 
shock model introduced in Sec. 3. They are more complete and realistic than first simulations of this model 
in [144]. Not only they include the effect of an external magnetic field, but more importantly they present 
the whole prompt emission during the time that a typical burst is observable. 

Each simulated burst consists of few regimes (at least 3): Beginning of the formation of an active region, the 
main emission, and a tail during which the emission becomes softer and fades off. Each regime corresponds 
to one of the models (3.26)-(3.29) for the evolution of width of the active region. All simulated bursts begin 
with steady state or exponential models because their initial width is zero. Some of the other parameters, 
specially indices parametrizing the evolution of parameters also can change with regime, but continuity 
conditions should be respected, see Sec 4.1. 

An important issue to be considered in numerical calculation of the synchrotron flux is the integration over 
Lorentz factor distribution of accelerated electrons in (3.32). The upper limit of this integral is oo, but in 
numerical calculations we must stop integration at a finite upper limit. It is easy to see that after a variable 
change the argument of the Bessel function in (3.32) becomes: 

y = J jx^J^ , LO'^ = Uj[\^,=^^ , X = -lehm ^ 1 (4.1) 

For 7e — )• oo, X — 7- oo and i^2/3(y) ~^ oo. Therefore, we must choose the upper limit of numerical integration 
such that at high energies where uj' /uj'^^ ^ 1, the variable x take much higher values such that y — )■ 0, 



otherwise the flux at high energies is under-estimated. For this reason we choose an energy dependent upper 
limit, Xmax = 200uj' / Lo!^. Evidently, when electrons distribution has an exponential cutoff, it automatically 
reduces the contribution of high 7e in the integral. Nonetheless, the above argument stays valid because the 
decreasing exponential and increasing Bessel function at ?/ — )■ cancel each others effect. 

4.1 Continuity conditions 

To have a smooth transition from one regime to another - the necessary condition for having continuous light 
curves - we have to respect some continuity conditions. Table 1 summarizes parameters of the simulated 
models. In each simulation all the relevant parameters should be given as initial values for the first regime. 
In the following regimes some of the parameters can be changed. They include the exponents/indices that 
determine the evolution of quantities such as and that respectively parametrize the evolution of the 
fractions of kinetic energy transferred to electric and magnetic fields Ce and cb, ctx the index determining 
variation with distance of the external magnetic field, Op index of power-law factor in electrons energy 
distribution, k index parametrizing the distance dependence of slow shell density, jcut the cutoff Lorentz 
factor of electrons when a power-law with exponential cutoff or broken power-law distributions is considered, 
and finally r, 6, and 5' that parametrize the evolution of active region width. As we explained before, the 
latter cannot be determined from first principles in this simplified model of relativistic shocks and needs 
input from detailed simulation of microphysics such as PIC simulations. 

Initial values of other quantities including the distance to the central engine rg, quantities that determine 
the kinematic of the shells such as 7'(rQ), r(ro) and their corresponding beta, the column density of the 
fast shell, the density of the slow shell, and the amplitude of the external magnetic field must be continu- 
ously transferred from one regime to the next. This means that their initial values in a new regime must 
correspond to their values at the end of the previous regime. Other quantities including the final T when 
shells are coalesced, the phase, and frequency of the external magnetic field stay the same in all regimes 
by definition. Note that although for the sake of simplicity we assume a constant precession period for the 
central engine and the external magnetic field, in reality they can change due to the slow-down of the engine 
and deformation of disk or envelop material that create the field. Nonetheless, if the slow-down time scale 
is much longer than burst duration this should be a valid assumption. If the engine collapses to a long life 
or a temporary magnetar before collapsing to a black hole, it is expected that its rotation rate and ejection 
of material continues roughly steadily for hundreds of seconds [135]. 

In addition to the matching of initial values, we must be also careful and use correct initial values for Tq, /3q, 
etc. which are important for the continuity of fiuxes. We should remind that evolution of physical quantities 
with distance/time in the model discussed here depend always on r' /t'q where Tq is the distance from the 
engine at the beginning of the regime - not at the beginning of the simulation. Therefore, in each regime 
we must use the initial value of r', /?', etc. of that regime in (3.22)-(3.24) and (3.26)-(3.29). 

4.2 Processes and parameters 

These simulations assume that the main energy dissipation process is synchrotron. According to previous 
analysis and simulations [19, 21] this is a reasonable assumption specially at middle range of energies. 
Compton scattering between high energy electrons and photons can be significant at very high or low 
energies, see next section, p — 7 scattering has been already shown to be insignificant and therefore its 
effect on the dynamics of the shock can be neglected [6, 18, 128, 129]. 7 — 7 — > e~e~^ can only affect 
photons of > 5 X lO^P ~ 5 X 10® eV. Moreover, most synchrotron photons are emitted in the direction 
perpendicular to the shock front and therefore are very close to collinear with each other. In this case the 
Center of Mass (CM) energy of photon pairs is very small, thus they cannot annihilate each other to e^. The 
situation is very different for a standard fireball model in which released photons have a high energy thermal 
distribution and their direction is random. As for the synchrotron self-absorption, it is only important at 
very low energies which are not crucial for the prompt emission. 
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Table 1 summarizes the parameters and initial values of physical quantities for the simulations presented in 
this section. Based on physical expectations not all of the parameters used in this model are independent. 
For instance, density, relative Lorentz factor of shells, and the fractions of kinetic energy transferred to 
electric and magnetic fields must be somehow related. However, the complexity and precision of plasma 
physics simulations do not yet allow a detailed exploration of parameter space to understand the underlying 
correlation between physical properties. Nonetheless, some of the relations between parameters are more 
evident. For instance, there is a degeneracy between the distance of the shock from the central engine 
and densities or column densities of shells that produce a given synchrotron flux, because due to adiabatic 
expansion the density of a shell decreases as even in absence of any energy dissipation. Therefore, 
larger the distance to the central engine, smaller the density of shells to generate the same flux. Other 
simple and obvious relations are: Emission hardness defined by a;^ is proportional to the magnetic field; 
Strength of the shock depends inversely on the quantity P defined in (3.25); Flux of synchrotron emission is 
proportional to the fraction of kinetic energy transferred to electrons; Duration of a burst depends on how 
quickly the induced magnetic field and accelerated electrons - equivalently electric field - approach to their 
maximum value and then decrease. We use these relations to design simulated bursts i.e. select initial value 
of parameters. Another guide in the selection of plausible values are observations. For instance, densities of 
shells are chosen based on the density of massive stars envelops, indices based on observations of GRBs, etc. 
Despite these relations, due to the large number of parameters in Table 1, a systematic exploration of the 
parameter space is not simple. Nonetheless, the consistency of simulations with observations significantly 
reduce the range of parameters. 

In the Introduction we raised the issue of radiation absorption inside the shocks. Here we address it in more 
details. High energy electrons are mostly ejected in front of the shock in slow shell [72, 111]. Assuming 
that the width of slow shell is roughly equal to the width of active region, the column density traversed by 
photons for the densest shells considered in these simulations is A'^^ ~ 10^^ cm~^ is ~ 10^^ cm~^/7' where 
7' is the Lorentz factor of the fast shell with respect to slow shell. The characteristic attenuation column 
density A for photons with energies in the range of 10 — 100 keV is A ~ 3 — 4 gr cm~^ or equivalently 
A ~ 2.5 X 10^^ cm~^ and larger at higher energies [73]. Therefore even for these high density shells, photons 
absorption inside the emission region is small < 25%. In the majority of our simulations the density of 
slow shell is much smaller and absorption is completely negligible. Evidently, the possibility of absorption 
by other material surrounding the central source cannot be ruled out. This issue is true irrespective of the 
distance to the central source. But, no strong absorption was detected in the GRB prompt emission. This 
is the evidence that surroundings of the source is cleared before GRB-making shocks occur - most probably 
by earlier ejecta. In fact the detection of precursors before the main peak and flares - softer emissions for 
thousands of seconds after gamma-ray spikes - is the evidence that activities of the central source last for 
much longer time than its peak activity observed as a GRB. Previous ejecta which were too soft to emit 
a gamma-ray flash most probably clear the surrounding of the source such that when the most energetic 
shells are ejected, there is no intercepting material around to attenuate them. Finally we remind that in 
equation (3.32) which describes the emission power, the value of r, Ar, r(r), and n'^ oc Nq are degenerate. 
Therefore, a larger distance to central source can be compensated by a smaller active region and lower 
densities for electrons. Thus, simulations presented here can be considered as representatives of a set of 
bursts with roughly similar x = '''o^^on'O^ f- The degeneracy is not exact because at longer distances the 
external magnetic field and thereby ujm are smaller. This leads to a softer emission. Nonetheless, it is always 
possible to adjust parameters to obtain a burst with the same characteristics occurring at larger distances. 
However, for longer distances from central engine, it would be more difficult - if not impossible - to explain 
the observed very short variabilities of emission. See also Sec. 4.7 for the issue of efficiency. 

Although at ejection jets are expected to be along the magnetic field lines, if the central source precesses at 
later times their direction can deviate from each other. The precession is due to internal forces and indepen- 
dent of the properties and behaviour of the jet after its release. Inside Alfven radius Ra ~ Rdisk acceleration 
and collimation is due to magneto-centrifugal forces that align the jet with the magnetic field [81]. Outside 
this zone kinetic energy of the jet dominates and collimation is most probably due to a helical induced 
magnetic field [1, 27]. In this case the change of magnetic field direction does not significantly change the 
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direction of ejected material, and in fact it helps its collimation by providing a transversal field compo- 
nent [1, 27]. Simulations show that progenitors of supernovae (collapsars) and their accretion disk indeed 
precess rapidly and irregularly [9] . Therefore it is expected that at the time of GRB formation a precessing 
magnetic field be present. Fig. 21 shows precession axis and its projection on the shock front surface. ^ 
is the angle between shock front and precession axis and p is the precession angle. They depend on the 
following quantities: direction and speed of the precession of the central source; precession of its magnetic 
field which can have an axis and period deviated from source ones; presence of a frozen field or a Poynting 
fiow [50]; distance of the shock front from the central source; etc.. The nature and value as well as relation 
between these quantities are mostly unknowns. For this reason we have considered arbitrary values for ^ 
and p. Also in all simulations their values are the same such that we can compare light curves without the 
infiuence of these poorly understood quantities. 

Apriori we do not have any rule for the duration of each regime. However, giving the fact that all evolving 
quantities depend on r/ro where rg is the initial radius of the current regime and the fact that /5 < 1, a 
given r/ro presents longer times for longer rg. This gives a rule of thumb for choosing the duration of a 
regime e.g. the initial rise of Eg and €b, etc. Another rule of thumb is the fact that if quantities become 
unphysically large, the solution of equation (3.22) leads to /3 < which is unphysical. This is an indicator 
that the value of evolving parameters are not anymore realistic. This seems very arbitrary, but sudden 
change of exponents in nonlinear systems most probably arise due to similar relations and internal feedback 
between quantities and conservation/ dynamical constraints that control their behaviours. Evidently, the 
ultimate constraint for us is obtaining simulated bursts which look like real ones, thus some try and error is 
necessary at this level. 

In [144] we have presented plots of various physical quantities involved in the model. For this reason we do 
not repeat them here and only present light curves and spectra of some examples of long and short simulated 
bursts which can be directly compared with observed GRBs. 

4.3 Simulated bursts: Light curves 

In this section we describe properties of light curves and spectra of simulated bursts. To simplify their 
comparison with each others and with real bursts we have transferred all long bursts to redshift z = 2 and 
all short bursts to z = 0.5 in a flat FLRW cosmology with Hq = 73 km sec~^ Mpc~^and = 0.73. 

4.3.1 Oscillations 

Between many simulations we have performed only a subset which are representative of the variety of 
observed bursts are discussed here. We remind that these simulations include only a single spike/peak of a 
burst. Well separated or even superimposed peaks in observed GRBs most probably originate from collision 
of shells which are separated in spacetime, or collision of successive shells with each others. In this case 
the accumulation of material in the slow-downed front shell leads to a stronger shock when lately emitted 
shells catch up and collide with it. Some of the brightest GRBs detected by the Swift Satellite such as GRB 
060105 [139], GRB 061007 [108, 109] and GRB 061121 [87] could be produced by such a process that we 
call highway collision because of its similarity to pile up collisions in highways. An analogue explanation 
is suggested for precursors that are seen in majority of bursts [29, 124]. On the other hand, even single or 
well separated peaks include substructures. The simplest explanation for their presence is inhomogeneities 
in the ejecta density. Nonetheless, our simulations show that substructures or even separate spikes can 
be produced in presence of a precessing magnetic field. In fact nonlinearity of the dynamics, dissipation 
processes, and fast variation of all physical quantities including the magnetic field can produce resonance 
emission, suppression, and/or create rapid variations in the emission i.e. spiky light curves with an acyclic 
appearance. A good example of such cases is simulation No. 1. in Table 1. Fig. 2 shows its simulated light 
curves without and with an external magnetic field. Without the presence of an external field as expected 
there is only one smooth peak in all energy channel. Moreover, lags between soft gamma-ray channels are 
somehow larger than what is usually observed in long bursts, i.e. up to few hundreds of milliseconds. When 
a precessing magnetic field with precession period of T = 2 sec (or equivalently / = 0.5 Hz) in the rest frame 
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Figure 2. Simulation No. 1 with no external field (lcft)and with \B\ = 35 kG (centre), \B\ = 70 kG (right). In all 
plots of light curves energy bands are written on the top of the plot and the color of their fonts correspond to the 
color of light curves of the band. 



of the engine is present, for the observer this burst includes 2 successive hard peaks and a third separated 
softer peak.^ There are also other peaks observable only in X-ray and softer channels. These peaks have 
different durations, and during about 12 sec that such a burst would be observable by a detector with 
sensitivity of the Swift-BAT, no periodic signature would be detectable in its PDS. Long observations in 
lower frequencies usually performed tens of seconds later i.e. at larger distances from central engine where 
the strength of the source field is significantly reduced. In fact, in this simulated burst there is practically no 
signature of oscillations in UV and optical because of the lag between their peaks. In addition, observations 
in these bands need at least few seconds of exposure and all data must be binned to remove/reduce the shot 
noise. This additionally smears the close to periodic behaviour of light curves. Note also that in presence of 
an external magnetic field the burst is harder, stronger and peaks have smaller lags. This burst looks like 
some of bright GRBs detected by Swift such as GRB 050525A and GRB 080409. In this simulation we have 
chosen relatively high densities to make a brightness comparable with very powerful GRBs such as GRB 
080319B and GRB 080607. Evidently weaker bursts are easier to simulate and a number of examples will 
be shown later. 

GRBs seem to have fast varying prompt gamma-ray radiation - a jitter - up to shortest revolvable time 
scales. It is very difficult to know if this high frequency component is intrinsic or due to the Poisson noise in 
the arrival time of photons. Nonetheless, various origins have been proposed for this phenomenon [63]. Here 
we argue that the origin of high frequency variations can be the precession of magnetic field of a forming 
magnetar or the magnetic field of the accretion disk of the engine during GRB emission. Fig. 3 shows 
the light curves of simulations with millisecond precession periods. It is evident that the precession of the 
field induces a fast varying component similar to the observations. Note that these light curves include 
substructures and overlapped peaks generated by variation of other quantities. Apriori both fast and slow 
varying magnetic fields can be simultaneously present. In fact it is well known that the differential rotation, 
turbulence, and instabilities in the accretion disk amplifies the seed magnetic field which in its turn produces 
a jet (see e.g. [50, 125]). This field can be the source for a slow varying component, coexisting with a 
fast varying component from e.g. the central source. There can be also a magnetic field frozen in the outer 
layers - convective layers - of the progenitor. In this case the total magnetic field that affects accelerated 
electrons in the shock is the sum of all these components. Fig. 4 shows an example of such situation. The 
external magnetic field energy is equally distributed in a slow component fi = 0.2 Hz and a fast component 
of /2 = 20 Hz. The signature of the fast component and its harmonics is clearly detectable in the PDS. Note 
that due to the dependence of ujm on the amplitude of the magnetic field, oscillation frequency is two times 
precession frequency. Only a non-significant peak at expected frequency of slow oscillation is detectable in 

'^For the sake of simplicity in all of the simulations presented here the precession frequency is the same in all the regimes. In 
reality we expect a variation due to acceleration of the central source after accretion of material or its slow-down due to energy 
loss by gravitational wave radiation or other forms of dissipation. 
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Figure 3. Top: Three examples of simulated burst with millisecond precession period at observe redshift. Simulation 
No. 2, / = 50 Hz, Tf = 100 (left), Simulation No. 3, / = 400 Hz (centre), Simulation No. 3, / = 50 Hz (right). 
Bottom: Ratio of external magnetic field amplitude to shock induced field for these simulated bursts. In these plots 
time is at the redshift of the source. 
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Figure 4. Light curves (left), differential power spectrum (center), and power spectrum of simulation No. 5 with 
a two component external magnetic field with rest frame frequencies of /i = 0.2 Hz and /2 = 20 Hz. Plots are at 
observer's redshift. Note that there is practically no signature of the long period magnetic field oscillations in PDS. 

the differential PDS. Fig. 4 shows that with available time resolution of detectors, even for bright bursts the 
high frequency part of PDS is dominated by Poisson noise and any signature of oscillation is smeared. Fig. 
3 includes also the evolution of ratio of external field to shock induced field for each simulated burst. From 
these plots it is evident that weaker the shock, more important the role of an external field in synchrotron 
emission. They also show how the field geometry affects its contribution at large distances and reduces the 
probability of detecting oscillations in the prompt low energy bands emission. We must also remind that 
in all the simulations presented here the strength of the external magnetic energy is orders of magnitude 
smaller than the kinetic energy of matter, thus bursts are not Poynting energy dominated. 

To be able to detect periodic features, the length of the interval during which the prompt emission is 
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Figure 5. Light curves (left), differential power spectrum (center), and power spectrum of simulation No. 2 with rest 
frame frequency /i = 0.5 Hz. 



observable must be long enough such that it covers a few cycles. An example of such bursts and their 
PDS is shown in Fig. 5. The precession period of the external magnetic field in this simulation is not very 
different from previous examples but the burst is longer - about 34 sec in gamma-ray channels. The light 
curves in all channels contain more features and a close to periodic behaviour is evident even by eye. But, 
it is noticeable that even in this case peaks are not completely periodic. Therefore, if the detection noise 
is high, a wide peak in the PDS of the burst can be interpreted as the effect of a white noise. In fact even 
in the noiseless PDS shown in Fig. 5 the peak due to the periodic component is not very significant. This 
explains the low significance of the peak in the PDS of GRB 090709A, if its origin is an external magnetic 
field. 

In conclusion, only under special conditions a periodic variation in the light curve of a GRB is observable. 
Notably, when the precession of the external magnetic field is very fast - a situation which is expected if 
the burst occurs after the formation of a temporary or long-life millisecond neutron star/magnetar - the 
signature of fast oscillations even for very bright bursts is strongly smeared by the Poisson noise and binning. 
Nonetheless, the clear detection of a periodic component in the light curve of GRB 090709a proves that 
maybe in this source the burst was formed before the rotation of the central source accelerates by accretion 
of material from the disk. Alternatively this relatively slow component can be due to the magnetic field 
of outer part of the disk or the envelop of the star, or formation of a long range density oscillation in the 
disk due to the presence of a magnetic field, see Appendix C for a brief review. As mentioned in Sec. 2.1, 
PDSs of many bursts seem to have weak periodic components, but due to the weakness of the signal and 
complexities described in the previous paragraph, these observations are not conclusive. In particular, as 
simulation examples here demonstrate, when burst duration is only few times the precession period, the 
effect of external magnetic field variation appears as substructures or simple peaks rather than coherent 
oscillations. This makes the detection of imprint of the precession very difficult. 

We single out two reasons for the absence of periodic features. The first one is intrinsic to the source. If 
oscillations have high frequencies, they would be smeared by the Poisson noise. Moreover, the geometry 
of the field plays a very important role both in the strength of the emission and in presence of oscillations 
specially in low energy bands. For instance, if the field is confined in a torus around the source, it affects 
the emission only at very early times when the jet traverses the torus. The second reason is the nonlinear 
dissipation and evolution of microphysics. Variation of indices that model the evolution of quantities such as 
€e, e_B , and densities along with nonlinearity of equations can produce sudden and resonant variations in the 
emission which breaks the coherence of slow precession of the magnetic field. Evidently, random variations 
of quantities due to turbulence in the jet environment also help to smear coherent oscillations. 

As we explained above, in the simulations presented here indices vary discretely. Although the variation of 
these quantities may be more gradual in the real bursts - equivalent to having a large number of regimes 
each lasting for a very short time - it is well known that in nonlinear systems phase transitions are usually 
first order i.e. discontinuous. Therefore, discontinuous variation of quantities implemented here can be an 
enough close imitation of real bursts. Moreover, a close inspection of simulated light curves show that they 
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do not produce the features explained above. In the next section we present more light curves in some of 
which the moment of regime change is visible and demonstrate that sudden changes in indices plays a minor 
role in creation of features. 

4.3.2 Phenomenological models 

In contrast to [144] simulations presented here cover the whole duration in which the prompt emission 
is observable by present instruments. Therefore, in this section we test the validity and differences of 
the phenomenological models (3.26)-(3.30) for the evolution of active region in internal shocks.® We also 
investigate the parameter space of the model, but due to its large volume we only present some examples 
and leave the systematic exploration and comparison with properties of candidate progenitors of gamma-ray 
burst to a future work. 

Fig. 6 shows light curves of 4 simulations. Each simulated burst includes 3 regimes. It is clear that when the 
middle regime is dynamical, the lags between gamma-ray bands are somehow large and peaks are softer. In 
absence of an external magnetic field irrespective of the middle regime the lags are too large to be consistent 
with observations. This proves that an external magnetic field can be crucial for generating bursts with 
observed properties. Nonetheless, considering the light curves in Fig. 2 this is not always the case and other 
physical properties of the shock are also important. This means that one cannot claim that the presence of 
an external magnetic field is a necessary condition for GRBs. In fact these simulations present very hard 
and bright bursts. Fig. 7 shows examples of less bright and softer bursts with and without the presence 
of an external magnetic field. In this figure simulations without an external magnetic field are consistent 
with observations. The case of the simulation in bottom-left plot of this figure is interesting. It is soft and 
consistent with late soft peaks detected in many bursts. Delayed bright peaks in X-ray bands can be detected 
as a bright X-ray flare without significant emission in gamma-ray. Many of this type of flares have been 
detected by the Swift-XRT instrument. In the bottom-right simulation rapid slope changes due to regime 
change are clearly visible and show that they are not responsible for the formation of features. We find that 
although the lags between various energy bands in these simulations are consistent with observations, they 



We leave the study of the propagation of internal shocks remnant to large distances and formation of external shocks with 
ISM/circumburst material to a future work. 
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Figure 7. Simulation No. 8 without magnetic field (top left) and with \B\ = 12 kG, / ~ 0.2 Hz (top right). Simulation 
No. 9 without magnetic field (bottom left) and with \B\ = 2.5 kG, / = 0.1 Hz (bottom right) 



are systematically larger. This means that at the beginning of the shock when it has its maximum strength, 
the synchrotron emission is hard and only when the shock is evolved and becomes softer, the emission in low 
energy bands becomes significant. This issue is directly related to the steeper spectral slope of simulated 
spectra in the framework of synchrotron model. In Sec. 4.4.1 we discuss this issue in details and relate it to 
the abrupt break in the distribution of Lorentz factor of electrons considered here. 

So far all the examples were long bursts. Short bursts are believed to be produced by the collision between 
two compact sources, e.g. 2 neutron stars, a neutron star and a black hole, etc. When one of the sources 
is a neutron star, one expects the presence of a strong magnetic field, specially at small distances from the 
central source. If the neutron star is young, the magnetic field on the surface of the star can be similar to 
magnetars \B\ > 10^^ Gauss. Even if the internal collisions occur at a distance of order 10^*^ cm and we 
assume that the magnetic field decreases as r~^, the external magnetic field would be few kGauss at the 
shock site and can have significant effect on the synchrotron emission. Fig. 8 shows the light curves of 4 
short simulated bursts. All these bursts have the same parameters except for the strength of the external 
magnetic field. It is evident that stronger fields create both brighter and harder bursts. Indeed without an 
external magnetic, according to the model presented here, this burst would be very weak and too soft to 
be consistent with observations. Lags also are shorter for stronger and harder bursts, although in the case 
of this example lags are very close to zero even in absence of an external magnetic field. By contrast, the 
total duration of the burst does not significantly depends on the strength of the external field. Because the 
precession period of the external magnetic field in this burst is only few times smaller than its duration, 
cyclic behaviour can be observed in the softest channels if simultaneous observations of high and low energy 
channels could be performed. To be able to better understand the effect of the precession speed of the 
external field in the features of short bursts Fig. 9 shows a simulated burst with the same parameters as 
Fig. 8 but a much faster precession / = 500 Hz. The comparison of light curves of these bursts show very 
little difference in the duration and strength of these burst, but when the precession is fast, the burst is 
slightly harder. We explain this observation by the fact that when the external magnetic field varies very 
quickly, the effect of having many field peaks at the initial part of the burst where it is harder leads to a 
slightly harder burst because a larger fraction of the available energy is emitted earlier, see the spectrum for 
each regime of this burst Fig. 14. In Fig. 8 we have also shown the ratio of the external to induced magnetic 
field. Given the fact that the distance of this burst from central engine as well as its duration are shorter 
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Figure 8. Simulation No. 10 without external magnetic field (top left), with \B\ ~ 10 kG, / ~ 5 Hz (top right), 
\B\ ^ 70 kG (bottom left), \B\ = 100 kG (bottom right). 




Figure 9. Light curves of simulation No. 10 with \B\ = 70 kG, / = 500 Hz (left); ratio of external to shock induced 
magnetic field (centre); light curves of the same model but with external magnetic evolving as (r/ro)^^ (right). 

than simulated bursts in Fig. 3, the external magnetic field is dominant in the whole duration of emission. 
Note also that in the middle of the burst when the shock is in its maximum strength, the contribution of 
external field decreases but it becomes dominant again close to the end when the shock becomes weaker. 
We also notice that the brightness and light curve of this simulation mimics closely the light curves of GRB 
05 1221 A the bright short GRB observed by the Swift-BAT. As Table 1 shows, simulations of this burst and 
its variants are performed by assuming a relatively flat magnetic field which varies only as (r/ro)~^. To see 
how much the geometry of the field can affect a short burst we have also simulated the same burst with 
{'i~/'''o)~^ during the whole burst. Fig. 9 shows also the light curves of this case. Their comparison shows 
that the burst becomes shorter and weaker. 

Fig. 10 shows another example of a simulated short burst with and without an external magnetic field. In 
contrast to the previous example this burst is enough hard even in absence of an external magnetic field 
to be consistent with observations. However, without an external field lags between the softest channels 
in gamma-ray are large. As we mentioned above, this must be related to energy distribution of electrons. 
The main difference between this burst and the previous one is the number density of protons and electrons 
in both slow and fast shells. They are 10 times higher than in the previous example. Consequently, the 
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Figure 10. Light curves of simulation No. 4 without an external magnetic field (left) and with \B\ = 35 kG, / — 500 
Hz (right). 



emission without an external field is stronger and harder. The precession frequency considered here is very 
high and consistent with the precession of a millisecond pulsar. We note that there are substructures in the 
light curves that do not exist in absence of the external field. As we mentioned before they are formed by 
nonlinearities, resonance effect, and rapid change of the magnetic field and evolution of indices. 

In these simulated short bursts the distance of shell collisions from central engine is at the shortest limit 
of the range given in the Introduction. Considering the large bulk Lorentz factor in these simulations, one 
obtains short variability time scales ~ At < 10~^ sec for them. A simple calculation shows that in a few 
other simulations presented here such as simulations No. 2 and 3 also At is very short because although they 
have larger distances, they have also large Lorentz factors. Nonetheless, they both generate long bursts, see 
Fig. 3. In particular, if simulation No. 2 occurs at 2; ~ 1, its light curve looks like a number of bursts with 
a duration of ~ 6 — 7 seconds which, depending on their redshift, can be considered as long or short. The 
lag in this simulation is also consistent with observations of latter category i.e. longer than short bursts 
but shorter than typical long bursts. In simulation No. 3 lags are longer and consistent with normal long 
bursts. Therefore, variability is not a determining parameter for classifying the progenitor of bursts to 
collision of compact sources and collapsar. As we ignore any relation between the radius of accretion disk in 
these sources and Lorentz factor of shells, the only constraint which can be imposed is similarity of emission 
to observed GRBs. At present, time resolution and sensitivity of detectors do not allow to distinguish 
between shot noise and genuine variability, thus any constraint on the latter is purely hypothetical. Even 
the constraint discussed in the Introduction is just an upper limit because inhomogeneities in the accretion 
disk can be very close and next to each others. We should also remind that as we discussed in Sec. 4.2, there 
is a degeneracy between the distance of the shocks from central engine and their density/column density. 
Thus, one can obtain a similar burst at larger distances i.e. with longer time variability, but its efficiency 
would be smaller, see Sec. 4.7 for more details. 

Considering the importance of an external magnetic field for synchrotron emission, the question arises 
whether a GRB can be formed without a shock. After all the internal shock model for GRB presumes shocks 
between shells that are accelerated to ultra-relativistic speeds by the engine through a mechanism which is 
not well understood but most probably is a combination of both mechanical and magnetohydrodynamical 
forces of an accretion disk. Due to the mass difference of electrons and protons it is expected that under the 
same force electrons are accelerated to much higher Lorentz factors. Therefore in presence of a transverse 
magnetic field they lose their energy by synchrotron emission. On the other hand, even in absence of a 
large scale precessing field a helical field is believed to be necessary for the collimation of jet outside Alfven 
radius where it does not anymore follow the parallel magnetic field [27] . This field dissipates electrons energy 
through synchrotron emission, and indeed recently this emission has been detected in numerical simulations 
of AGNs jets [98]. The emission is expected to be low energy - in radio band for AGNs jets with F '-^ 10 
and in IR/optical band for ultra-relativistic jets of GRBs with F ~ 1000, assuming the same strength for 
the helical field. Such radiation from GRBs has never been observed. The reason is most probably the 
absorption of this emission be the jet itself because such emission occurs along the jet. Relatively low 
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Figure 11. Light curves of simulation No. 11 witliout an external magnetic field (left) and with \B\ — 100 G, f — 0.2 
Hz (right). 



energy photons are absorbed or scattered during their propagation through the large column density of the 
jet. Moreover, electrons lose their energy very quickly and in a flash. In a collisionless flow they are not 
reaccelerated, and are simply dragged forward by protons flow to which they are electrically coupled. In this 
case if at later times the external magnetic field changes its direction no accelerated electron is available. 
Only when electrons are reaccelerated in the shock, the external transversal field can become active again 
and dissipate electrons energy by synchrotron emission. Shocks accelerate electrons are to higher Lorentz 
factors and generate additional magnetic field, thus the synchrotron emission is in much higher energy bands. 
Moreover, a shock can sustain a continuous acceleration of electrons for a relatively long time by gradually 
transferring the kinetic energy of colliding baryons to electrons. Therefore, even in the situation in which 
the magnetic field made by the shock is subdominant, the presence of a shock is a necessary condition for 
the formation of a GRB. In Poynting flow models reconnection of magnetic lines plays the role of shocks for 
accelerating electron. But, as mentioned before, the long distance of the process from the engine and rarity 
of line connection disfavour this model. 



Finally we notice that in all the simulations presented here the value of eg 



< 0.15. This is consistent 



with the values obtained in numerical simulations of electron acceleration in PIC [111] and Monte Carlo 
simulations [24] (see [8] for a review). On the other hand, this value is orders of magnitude larger than 
^e\max ~ 10~^ found in some modification of basic SSC model to explain the behaviour of high energy 
spectrum of GRBs [121]. In addition to a low efficiency, this type of models produce a significant bump at 
~ 1 GeV which has not been observed. We should also remind that eg, density of shells, and distance 
to the central source are somehow degenerate i.e. a burst with roughly similar properties can be obtained 
for eg ~ 0.3, rg ~ 10^^ cm and smaller densities. Fig. 11 shows an example of such case for simulation No. 
11. The reason for such degeneracies is clear. At large distances densities are smaller and therefore higher 
efficiency is necessary to produce the same amount of photons. At smaller distances a large Ce leads to too 
efficient and hard emission and the duration of burst is usually very short, only consistent with very short 
hard bursts. Although apriori large eg can be consistent with observations, giving smaller densities at larger 
distances, they are unrealistic. Moreover, the peak in the light curves shown in Fig. 11 is too long to be 
consistent with a usual GRB, although GRBs for which a supernova is observed e.g. GRB 060218 [110] are 
very faint and long. In contrast to other GRBs, they usually have a thermal spectrum with a temperature 
of few keV. The gamma-ray is in fact the high energy wing of this thermal emission which most probably 
originates from an uncollimated - or mildly collimated - ejecta. The lack of significant variability in such 
bursts and simulations here is another argument for occurrence of normal GRBs at much shorter distances. 



4.3.3 Tail emission 

In real GRBs usually a long emission tail specially in X-ray is observed at the end of the main spikes 
with a very steep power-law or in some cases exponential decline after few hundreds of seconds. In these 
simulations we see the steep decline, see Fig. 17 where light curves are plotted in logarithmic scales. 
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Figure 12. Light curves (left) and spectra (riglit) of simulation No. 12 which simulates weak shocks or when the 
main shock is weakened due to close to complete coalescence of two shells. In right plot the highest full line curve is 
the total spectrum and other curves are spectra in the time interval of 4 simulated regimes. They show how the peak 
of spectrum moves toward lower energies when the shock becomes weaker. 



This figure also shows that the decline slope is flatter in lower energy bands which is consistent with the 
observations, and explains why X-ray and optical light curves of GRBs have different slopes and breaks 
are chromatic. However, in our simulations the steep decline occurs much earlier than in real GRBs. We 
have two explanations for this difference. In what concerns simulations, we consider only one shock with 
few regimes i.e. we end the shock quite quickly. In fact due to numerical complexities and validity range 
of approximations when the shock becomes very weak, simulations are usually stopped before /3' ~ i.e. 
before shells coalesce completely. In real bursts one expects that shocks fades more slowly, thus the soft 
emission continues for much longer time. The second reason is intrinsic to the activity of the engine. Usually 
late peaks are softer meaning that gradually the process responsible for the acceleration of the jet becomes 
weaker and less efficient. Therefore the tail of the prompt emission is due to soft shocks. 

Rather than following slow decay of simulated prompt emissions, we simulate an example of weak shocks 
using simulation No. 12 which is similar to simulation No. 11 but has a longer distance from central source 
and lower densities. The relative Lorentz factor is also much smaller to simulate weak shocks when shells 
are very close to complete their coalescence. Light curves and spectrum of this simulation are shown in Fig. 
12. The first few seconds when shocks is formed must be ignored because we had to make the shock from 
beginning. The tail emission presents the declining regime of the light curves. We notice that similar to 
what is observed in many bursts, they are smooth and have a steep break after few hundreds of seconds. 
They confirm the conclusion taken from simulation No. 11 that longer the distance of the shock from the 
central source, smoother light curves and softer the burst. We also remind that these simulations do not 
take into account the curvature of the emission surface. It somehow flattens the break slope due to later 
arrival time of photons. In radiative shocks the steep break occurs when the two shells are coalesced. In 
the example presented here, the simulation is stopped when the Lorentz factor of the fast shell reduces to 
400.002 for an input value for the final Lorentz factor of 400, see Table 1. 

This simulation is consistent with the results of [131] who relate the steep decline and the plateau regime 
of the X-ray flux respectively to the tail of prompt emission and its overlap with the rise of the afterglow 
(which is not simulated here). The peak of emission spectrum for this shock is at few hundreds of eV, much 
lower than in typical prompt emissions. 



4.4 Simulated bursts: Spectra 



As we mentioned in the Introduction, since the launch of Fermi satellite and observation of prompt emission 
of GRBs in MeV and GeV energy bands, various models are suggested to explain unexpected results of 
these observations. Thus, an evident validity test for the model presented in [143] and extended here is its 
ability to reproduce the observed broad band spectra, specially at high energies. In the context of internal 



shock SSC model, [13] and [22] have aheady performed detailed simulations of GRBs spectra based on a 
formulation of shocks by [19]. They find a bump in the spectrum due to the inverse Compton scattering 
between high energy electrons and photons at E lOGeV, following by an exponential break. According 
to their predictions and depending on the parameters of their model, other features can also be present in 
the spectra at high energies. None of these features have been observed in the Fermi-L AT/GBM data. By 
contrast, in some bursts such as GRB 090510 it seems that the high energy part of the spectrum up to 
highest observed energies is best fitted by a power-law with positive index. 

In this section include only present the spectrum of synchrotron emission in simulated bursts. Compton 
emission is discussed separately in Sec. 4.5 and we find that in most cases it has a negligible contribution to 
the total emission in GeV energy band. Because the spectra of simulated bursts discussed in the previous 
section have very similar properties we do not present the spectra of all the simulations and concentrate 
only on some of them which are representative. 

Fig. 13 shows the spectra of simulation No. 7 during separately for each regime as well as the total spectrum 
in observer frame. Each plot includes spectra for a set of index and exponential cutoff values. The peak 
of photon number per unit energy N{E) appears as a bump with an amplitude roughly proportional to 
the index p of the power-law factor of electrons Lorentz factor distribution.^ It is directly related to the 
average synchrotron characteristic frequency in the time interval in which the spectrum is determined [47] 
and its position determines the hardness of the burst. When Ucut/^m ~ 10 — 100 the peak is followed by 
another bump at higher energies which is due to the exponential cutoff in electrons energy distribution. 
A good example of such bursts is GRB 090926A [33]. When p < 2 and cocut/^m ^ few 100, the bump 
becomes flatten and looks like a linear increase in the emission similar to behaviour of the spectrum of 
GRB 090510 [84, 96]. These simulations show that this feature is due to the high value of cutoff in the 
electrons energy distribution and no additional process is needed to explain it. The spectrum of this burst 
also shows clearly the low energy peak around ~ MeV energies. When the cutoff is low or p > 2, this peak 
becomes the absolute peak of the fluence. As Fig. 13 clearly shows its position evolves with time. This 
is consistent with the observed evolution of spectral peak GRBs observed by the Fermi-LAT such as GRB 
080916C [30, 43, 114]. 

in these simulation at E <^ Ep^ak the slope of the spectra approach the theoretical value of 4/3. However, 
in the observed GRBs low energy slope of gamma-ray emission is usually smaller than this value and is on 
average ~ 1 [104, 138]. Addition of inverse Compton emission slightly flatten the spectrum both at low and 
high wing of this simulated burst, see Fig. 17. This is in agreement with simulations of [22]. However, in 
contrast to some of previous simulations we find that in very hard bursts such as GRB 090510 Compton 
scattering does not have a significant contribution in the high energy emission. This explains the lack of 
any feature due to this process in the spectrum of GRB 090510 which was a short hard burst. See the 
next section for more discussion about low energy spectrum. Although the high energy spectra of some 
bursts are consistent with an exponential cutoff, many bursts show broader spectra. Such a behaviour can 
be reproduce be considering a broken power-law distribution for electrons at high energies rather than an 
exponential cutoff. Few examples such bursts are shown in Fig 15-a & b. It is clear that in this case the 
spectrum is more featureless and broader. Finally we notice some features in the low energy wings of the 
total spectrum of some simulated bursts - see Fig. 14. Our tests show that these features are artifacts of 
the abrupt change of parameters from one regime to another. Fig. 15-c shows two simulations with same 
parameters but one of them (dash line) has a larger number of regimes and variation of parameters from 
one regime to the next is smaller. The spectrum of this simulation does not show any low-energy feature. 

We should also remind that in all the simulations presented here the selection of electron index and cutoff is 
completely arbitrary and no connection between these parameters and other quantities that determine the 
physics of the jet and shock are considered, simply because we ignore any relation. Some PIC simulations 
confirm a power-law with exponential cutoff distribution for electrons [111], but do not give any hint about 

^This bump is not exactly at the same energy as the peak of N{E) but they are close and from now on we call it the peak 
although it is not always the absolute peak of the fluence E'^N{E) which are plotted here. The reason for plotting the latter 
rather than N{E) is that spectral features are enhanced in E^N{E). 
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Figure 13. Simulation No. 7: a), b) and c) are respectively spectra of the first, second and third regimes, d) is the 
total spectrum. In each plot curves correspond to different values for p and LOcut/'^m'- P = 2.5 and ujcut/'^rn — 0.5 (full 
line), 1 (dash-3 dots), 10 (dot), 100 (dash-dot); p = 1.9 and ijjcut/'^m = 1000 (dash). The external magnetic field in 
these simulations is 10 kG. The low amplitude full line has p = 1.9 and ujcut/^m ~ 1000 but no external magnetic 
field. 
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Figure 14. Simulation No. 10 with power-law -I- exponential cutoff for electrons, from left to right spectra of 3 
regimes and total spectrum. Top row: Simulations with various values of power-law index and cutoff and |Bea:t| = 70 
KG: Ucut/^m = 1000, p = 1.5 (full line); uJcut/^m = 100, p = 1.5 (dash); uicut/^m = 3, p = 2 (dot-dash); uscut/^^m = 3, 
p = 2.5 (dot). Bottom row: Simulations with various values of parameters and iSextl = 100: Wcut/wm = 1000, p ~ 2.5, 
Eg — 0.002, r — 500 (full line); utcut/^m and p as previous case and ~ 0.02, F = 50 (dash); Uq ~ 5 x 10^^ cm"'^ and 
other parameters as the previous case (dot-dash); varying p with index —0.2, 0, 0.5, initial p = 2.5 and Ce = 0.002 
(dot); initial p = 1.8, the same indices as previous, and uJcut/^m = 0.5, 1000, 100 (dash-3 dots). Different spectral 
shape of the last simulation in different time intervals is similar to observed spectra of GRB 090510. Parameters which 
are not specified here are equal to what is given in Table 1. 
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Figure 15. a) Simulation No. 7 with a broken slope at ujcutl'-^m = 3, pi = 2.1, P2 = 4 (full line); p = 2.5 (dash 
line); b) Simulation No. 10 with a broken slope at lOcutl'-^m a) and pi = 2.5, ■p2 = 4, j-Bertl = 17 kGauss (full line); 
Pi = 2.1, p2 = 4 and same |i3ea;t| as previous (dash); pi = 2.1, -pi = 3, |i?ea:t| = 26 kGauss (dot-dash); same slope 
and \Bf.xt\ — 35 kGauss (dot); same slope and |Be£ct| = 70 kGauss (dash-3 dots), c) Simulation No. 5 with same 
parameters: 4 regimes (full line), 7 regimes (dash). With more gradual change of parameters in the latter case the 
spectrum at low energies becomes smooth and the low energy feature disappears. 



the relation between the value of index and cutoff energy and other quantities. Considering the random 
nature of instabilities, it is plausible that feedback between various quantities and the randomness of the 
latter smear any feature and create the observed featureless spectra - see also below for more explanation. 

4.4.1 spectrum slope at Xgsn energies 

As we mentioned above, the spectrum of simulated bursts at energies much smaller than peak energy - at 
least two or three orders of magnitude - is consistent with the theoretical prediction of ~ 4/3. This is a 
well known result. From properties of synchrotron emission [97, 99, 105, 106], in another word the Bessel 
functions in equation (3.32), this slope is obtained for uj' <C Wc- However, the average low energy slope 
of observed long bursts is about 1 [104], although brighter bursts are somehow harder and have steeper 
spectra at low energies [48]. Such a flat slope can be explained if electrons are in what is called slow cooling 
regime [106] (see also Appendix D for more details), but the efficiency of synchrotron emission in this regime 
is low. For this reason the general believe is that in GRBs electrons do not cool slowly and are in fast cooling 
regime. Thus only steep slopes close to the theoretical value of 4/3 are expected if the origin of emission 
is synchrotron. A number of solutions such as the effect of Klein-Nishina and adiabatic expansion [22] as 
well as thermal emission have been suggested to explain the flatness of low energy spectrum. These works 
conclude that synchrotron self-Compton processes can explain observations although for GRBs with very 
flat low energy spectra another origin, e.g. thermal emission may be necessary. 

In this section we consider this issue in the framework of the formulation of synchrotron discussed in [143] and 
simulated in this work. On the theoretical side we should remind that the strict classification of accelerated 
electrons in a shock to slow or fast cooling by comparing 7^ with cooling Lorentz factor 7c defined in (D.l) 
is meaningful only for a power-law distribution with an abrupt break at 7^. In this case most of electrons 
have a Lorentz factor close to 7m, thus comparison of 7^ with 7c indicates whether majority of photons are 
emitted with energies below the peak or above it. If the distribution of 7e is such that there is no abrupt 
break but it has a tail of low energy electrons that preferentially emit synchrotron photons at low energies, 
a flat spectrum in both sides of the peak energy can be obtained. Because our theoretical calculation is 
performed only for simple power-law or power-law with exponential break, here we try to simulate more 
complex simulations by considering a variable distribution for electrons in different regimes of a burst. 

Spectra presented in Figs. 12, 13, 14, 15 show that with evolution of a burst, the peak energy gradually 
moves to lower energies. The spectrum at high energies can have variety of behaviour, notably it can be 
very flat or even rising if electrons spectrum is shallow. Fig. 16 shows the spectra of the simulation No. 1 
with various spectral behaviour for electrons. The slope of the fluence in these simulations is < 1 from the 
peak up to ~ 1 keV for Fig.l6-a and up to ~ 10 keV for Fig.l6-b which has much larger external magnetic 
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Figure 16. Spectra of simulation No. 1 at £^ < Epeak in 4 regimes/time intervals with power-law-exponential cutoff 
distribution for electrons in the first and last regimes and broken power-law in two middle regimes. The full line is 
the total spectrum, a) Electron distribution 4 regimes: 1- power-law with exponential cutoff p = 0.5, ^cut/lm — 10; 
2- broken power-law pi = 1.5, p2 = 2.5, ^cutllm = 10; 3- broken power-law pi = 2.5, p2 = 3, ^cut/lm = 10; 4- 
power-law with exponential cutoff p = 1, ^cut/lm = 10; |5ea;t| = 10 kG. b) 1- power-law with exponential cutoff 
p ~ 0.5, ^cut/lm = 10; 2- broken power-law pi = 1.5, p2 = 2.2, ^cut/lm — 5; 3- broken power-law pi ~ 2.2, p2 = 
5, "fcut/jm = 10; 4- power-law with exponential cutoff p = 1, "fcut/"fm = 1; |^ea;t| = 100. Note that in both cases the 
slope of the spectra up to ^ 2 orders of magnitude from the peak are flatter than 4/3 and break to this value at very 
low energies. 



field and is harder. Luminosity of the simulation in Fig.l6-a is not very strong. Therefore such a model can 
not represent a typical GRB, but the emission from a sub-population of mildly accelerated electrons which 
are present at the same time as highly accelerated electron making the high energy component. Simulation 
in Fig.l6-b is hard and presents the spectrum of a typical GRB. 

In exceptional cases such as in GRB 061121 [87] it was possible to observe simultaneously the burst in 
gamma-ray and x-ray. Similar to these simulated bursts, the slope at low energy wing of GRB 061121 is 
shallow and approaches 4/3 only in soft x-ray and optical/IR [88]. 

4.4.2 Electrons distribution 

In the previous sub-section we showed that the presence of mildly accelerated electrons can produce a flat 
low energy spectrum. Does such a population of low energy electrons exist ? Both PIC [111] and Monte 
Carlo simulations [24] of charged particle acceleration in relativistic shocks show that only about 10% of 
electrons are accelerated to very high energies. The rest of electrons keep their original thermal distribu- 
tion. But no signature of a thermal emission is found in the observations. In presence of the magnetic field 
induced by the shock they emit synchrotron radiation which can contribute to low energy emission and 
notably flatten the spectrum. On the other hand, it seems unlikely that such a large fraction of electrons 
stay unmodified. Moreover, the absence of any thermal signature means that either even before the shock 
they had a nonthermal distribution, or their distribution somehow has been modified during the shock. In 
fact spectra of GRB 080916C which was enough bright to allow the determination of the broad band spectra 
in different time intervals [43] show that the low energy slope of the first interval between trigger until about 
7 sec is much steeper than later intervals [114]. We interpret this as the presence of an abrupt break, similar 
to simulations here, in the electron distribution during the first time interval. By contrast, in the second 
interval the break is smoother meaning that a larger number of electrons are moderately accelerated and 
their distribution at low energies is flatter. Considering the results of PIC simulations of [72], the micro- 
physics of these observations can be qualitatively explained as the following: 

During the first interval the amplitude of the electromagnetic wave is not yet very high and free path length 
of electrons scattered by their local environment i.e. protons/ions and other electrons \ < Vg = m'yf./eB 
gyro-radius (Lamb radius). Therefore, before these electrons can diffuse and accelerate, they are scattered 

by other particles. It is why simulations find them in the same distribution as protons which are practically 
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unperturbed by the collision [24, 72]. Gradually the amplitude of electromagnetic wave increases, thus more 
electrons are accelerated and lose their energy by synchrotron emission before being scattered by their envi- 
ronment. These mildly accelerated electrons which necessarily emit low energy photons play an important 
role in smoothing electrons Lorentz factor distribution, and thereby flattening the photon distribution at 
low energies. This should also make smaller lag between high energy and low energy band, consistent with 
observations. In conclusion, we should use smoother distributions at low energies for electrons. 



4.5 Compton scattering 

Inverse Compton (IC) emission from GRBs is studied by various authors see e.g. [13]. In this section we 
briefly describe Compton scattering in the framework of model and simulations presented here. As the 
interaction between photons and accelerated electrons does not significantly affect dynamics of the shock 
we can determine the flux of IC and its spectrum separately. We should also remind that for far observers 
both electrons and synchrotron photons are boosted in their direction. Therefore, the Center of Mass (CM) 
energy of their interactions is small and very little energy is exchanged between them. This qualitative 
argument shows that inverse Compton scattering cannot be a significant source of high energy photons. 
By contrast, due to delay between production of photons by synchrotron emission and their scattering by 
electrons, their light curves are more extended in time. Simulations discussed below are consistent with this 
description. 

Using the Boltzmann equation, the rate of Compton scattering can be written as [141]: 



dt' J^^ 7e J E[ ^' ^'^'^'J ' ' ' E'dE'de' 

A{je,E[)^Pi.P; (4.2) 

where F'^^^ and F'^^ = dPgy/uj'duj' are respectively the number density of photons produced by Compton 
scattering and photons produced by synchrotron emission, 9' is the angle between the scattered photon and 
its original direction, and are respectively 4-momentum vectors of incoming electron and photon. 
We have assumed head-to-head collision to obtain an upper limit on the Compton scattered photons and 
integrated over all possible directions in the rest frame of the shock front/active region. Both Thompson 
and Klein-Nishina terms (s and t channels) are included in the Compton cross-section a. To determine 
F'^\ one has to integrate the right hand side of (4.2) with respect to time. However, as the final aim is to 
calculate the flux seen by a far observer, before integrating with respect to 6' one has to multiply the right 
hand side by Tf{l + /3/cos(0') [103] to take into account angular transformation with respect to observer. 
Another important fact to be considered is the extent of distribution of accelerated electrons inside the slow 
shell. As explained in Sec. ??, synchrotron photons are emitted by accelerated electrons in the Energetic 
Electromagnetic Structure (EES) very close to shock front. PIC simulations show that the width of EES is 
only ~ 10 times of Agp the plasma wavelength of electrons, and accelerated electrons diffuse inside slow shell 
~ lOOOAep [72, 111]. Therefore synchrotron photons that are not scattered in this thin layer, do not meet any 
accelerated electron coming toward them (in the slow shell frame). For this reason the integration of (4.2) 
must be limited to the time that photons pass through this region. If one naively integrates this equation 
for the full passage time through slow shell - assuming that its end coincides with the end of the synchrotron 
emission - one obtains, in contradiction with observations [43], a significant deformation of spectrum from 
synchrotron. Compton parameter for both cases are shown in Fig. 17-e. Evidently, the absence of concrete 
information about micro-physics of the shock create a large uncertainty on the determination of Compton 
parameter and thereby, the contribution of Compton scattering in shaping the observed spectrum. Therefore 
simulations reported here must be considered as a rough estimation. 

For angular cross-section we have used full relativistic quantum expression [56]. Fig. 17 shows the light 
curves, spectrum, and Compton parameter of Compton scattered photons in simulation No. 1. For com- 
parison the light curves and total spectrum of synchrotron alone are also shown. As expected, Compton 

emission decreases with time much slower than synchrotron emission. Nonetheless, it stays - specially at 

- 30 - 



high energies - orders of magnitude smaller than synchrotron flux. This would not be the case if we as- 
sume a uniform distribution of accelerated electrons everywhere along synchrotron outgoing path, see large 
Compton parameter for this case in Fig. 17-e. Therefore, according to these simulations inclusion of inverse 
Compton scattering cannot explain the high energy tail of emission observed by the Fermi-LAT. These re- 
sults are consistent with the conclusion of [43] that finds only one process responsible for the emission in all 
energy ranges. The suggestion of [22] that Klein-Nishina inverse Compton cooling can explain the shallow 
slope of Fermi and Swift observations at E < Ep(.ak can be true only if the penetration length of accelerated 
electrons inside slow shell is much longer than what PIC simulations obtain. However, this will also change 
the spectral shape at low energies and should produce a bright emission in UV/optical bands unless these 
photons are absorbed. Nonetheless, a small contribution from Compton scattering can somehow flatten low 
energy spectrum and makes it more consistent with observations. As we suggested in Sec. 4.4.1, a flat 
electron distribution is most probably necessary to explain the high energy slope. In the next section we 
suggest an alternative explanation for the observed high energy tail emission. 

4.6 Origin of the high energy tail 

For having a strong Compton emission, accelerated electrons moving toward the observer must collide - in 
preference head-to-head - with photons emitted by a source at rest or slow moving, for instance photons 
from material surrounding the source. The problem with this suggestion is that at late times the shock is 
much weaker, thus not sufficient high energy electrons are around. Moreover, at these late times the jet 
is well outside the star envelop, thus the density of environment photons is relatively low. A number of 
other explanations are suggested for the tail emission. For instance, external forward shocks [25, 28, 37, 38, 
54, 55, 70, 71, 122, 145], Compton scattering of photospheric photons by upstream scattered high energy 
electrons [123], etc. Although both of these processes can have some contribution in the observed flux, it 
seems unlikely that they provide a satisfactory explanation for all the high energy bursts observed by the 
Fermi-LAT. External shocks are too weak to explain the observed flux [137] and its fast variabilities [123]. 
Interaction of up-scattered electrons is evidently a possibility, but it is originally suggested in the context of 
standard fireball models with photospheric emission before shocks. It needs some fine-tuning on the delay 
between two successive shells (or fireballs) and distance between shocks and photosphere. Moreover, the 
peak of high energy emission must be simultaneous with a later peak at lower energies, but the peak of high 
energy light curve does not always correspond to a peak at lower energies, see for instance GRB 090902B [32] 
where 2 late peaks in the Fermi-LAT light curve of photons with E > 100 MeV do not correspond to any 
peak at lower energies. 

Comparing simulated light curves and spectra in Fig. 17 with observations, we find that synchrotron 
emission can explain the main high energy emission which in all the observed bursts coincides with early 
peaks in lower energies. On the other hand, it seems that in bursts with long high energy tails later peaks 
are significantly harder than earlier ones. Either they have occurred at the same time as earlier ones but 
at shorter distances from the engine, or more probably they are due to the collision of late ejected shells 
with the remnant of the earlier collisions which are slowed down. In this case their hardness can be due to 
the accumulation of material in the slowed down earlier shells, what we called highway effect in Sec. 4.3.1. 
As we showed in the previous section, although inverse Compton scattering last for a much longer time, its 
flux is not enough to explain the tail emission. If there are additional photons from another source rather 
than synchrotron emission, the efficiency of Compton scattering would increase. However, no signature of 
these photons which must also contribute to low energy emission has been found but it is possible that they 
can be absorbed [95]. Indeed, only very few observations of the prompt emission in low energy bands have 
been performed and the quality of data does not allow to estimate the column density of material and the 
absorption at the time of prompt emission. An exceptional case is GRB 061121 [87] which its main peak 
occurred about 70 sec after the precursor. Its average column density during the peak is estimated to be 
Nh ~ 1.4 X 10^^ but Swift-XRT data show fast evolution [88], and therefore the initial column density can 
be much larger. In any case, it does not seem that this column density be enough to completely absorb low 
energy emission and in fact the afterglow of this burst was observed in optical and IR [87] . 

Here we suggest that high energy tail emission can be due to a delayed synchrotron emission of the most 
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Figure 17. Synchrotron light curves (a) and spectrum (b), inverse Compton light curves (c) and spectrum (d) 
for simulation No. 1 with p ~ 2.5 and ^cut/lm = 1000. Compton parameter (e): assuming the same density of 
accelerated electron everywhere in the slow shell (blue), for the layer containing accelerated electrons according to 
PIC simulations (magenta). We have neglected scattering/absorption of photons by cold plasma (see Sec. 4.2) in 
calculation of Compton parameter. Compton scattering light curves and spectrum correspond to the latter case. 
Contribution of t-channel (f ) with respect to energy band of the scattered photon. Each x symbol indicates the ratio 
of t-channel to total in the right hand side of equation (4.2) at a given time. Sudden break at low energies is due to 
truncated synchrotron spectrum in our simulations. 

energetic electrons. To be accelerated, electrons must oscillate back and forth in the collided region - the 

active region - where the electromagnetic wave is concentrated [24, 72, 111]. However, due to the strength 

of normal field their movement is oblique to the flow direction and helical, and the minority of them which 

are accelerated to very high energies have also a very large Larmor radius. With the progress of the shock 

gradually the electric and magnetic fields inside the active region would not be enough to confine these high 

energy particles and they escape to the down stream where the magnetic field flux is very small, see Fig 1 of 

[111]. A similar process is believed to be the origin of cosmic rays in the shocks at much larger spatial scales 

- see e.g. [3]. These particles do not lose their energy until either are Compton scattered by photons or in 

presence of an external magnetic field dissipate their energy as synchrotron photon. We must also remind 
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that the highly charged region in the shock front [111] at least partially screens the down stream from the 
magnetic field of the source. This process has been studied in the context of accretion from a disk in which 
the accretion flow partially screens the disk from the magnetic field of the central source, see e.g. [57]. Such 
a reservoir for high energy electrons and their delayed dissipation explains why the spectrum of GRBs is 
consistent with one component emission [43] but their light curves seem to show two components. 

To estimate the delay time we assume that the energy of observed high energy photons is close to the syn- 
chrotron characteristic frequency of the emitting electrons. Using (3.31) and the characteristic synchrotron 
emission time of electrons r' = cme7e/ei?' = 3tt^^/uj'^, and a Lorentz factor of F ~ 1000, we find B' < 300 
Gauss to have a maximum delay of ~ 100 sec for photons with E ~ lO^'' eV for the observer. This means 
that the environment around the shock front must be either very unmagnetized or very efficiently screened 
from the magnetic field. 

Another possibility is a resonant movement of the most energetic electrons with Electromagnetic Energy 
Structure (EES) that forms in the two sides of the shock front [72, 111]. The maximum of electrostatic and 
electromagnetic fields are situated very close to the shock front and electrons are accelerated by moving back 
and forth between down stream and up stream in this field. The electromagnetic component is dominated 
by transverse electric and magnetic fields and is approximately a standing wave expanding more in the 
downstream than upstream, see figure 2 of [72] and figure 3 of [111]. Due to the phase difference between 
electric and magnetic fields, electrons are accelerated where the electric field is strong and lose energy where 
the magnetic field is strong. Although we do not yet have a complete simulation of this process, it is 
perceivable that with the evolution of the shock the electromagnetic field gradually becomes non-standing 
and propagates to down stream. Electrons with highest Lorentz factors can follow the propagation of 
the field for sometime and ones which are trapped in the electric field dominated region continue to be 
accelerated until their lag from the electromagnetic field bring them to the magnetic field dominated region 
where they lose their energy by synchrotron emission. This phenomenon also explains why a tail emission 
in low energies is not observed. Low accelerated electrons which emit in lower energies cannot follow the 
propagation of the electromagnetic wave (EES) and their lag with respect to EES take them quickly to 
the the magnetic field dominated region where they dissipate their energy quickly. Figure 2 of [72] shows 
that the wavelength of EES A'^;^;^ ~ SA'^^^^ where A'^^^^ = c/Up and uj'p = {en' /mpeoY^"^ are respectively 
skin depth and plasma characteristic frequency of protons in the fast shell. Using the same typical values 
for these quantities as what we used in the previous paragraph, n' ~ 10^^ cm~^, and assuming that EES 
propagate with a speed very close to light propagation in vacuum, we find a time delay of ~ 800 sec for 
£^ > 10 photons. Evidently, this is just an order of magnitude estimation because PIC simulations used here 
are not realistic for a GRB and we ignore the exact evolution of EES. Nonetheless, this exercises shows that 
in principle this idea is viable, and more realistic PIC simulations should be able to verify it. 

4.7 Efficiency 

The efficiency of synchrotron self-Compton models has been a serious issue in the internal shock mod- 
els [137]. PIC simulations find that only ~ 10% of the available kinetic energy is transferred to accelerated 
electrons [111]. Table 1 shows that in our phenomenological simulations also eg is always small and its 
maximum value is always ^ 10 — 15%. For higher values the duration of the burst becomes very short. 
The second issue is how much of the energy transferred to electrons is effectively emitted as high energy 
synchrotron. In absence of an external magnetic field the efficiency of emission depends on e^. PIC simula- 
tions find < 0.001 for a baryonic plasma. Although higher values apriori can increase the efficiency, they 
also shift the peak of spectrum to very high energies inconsistent with observations, see equations (3.17) 
and (3.31). According to (3.5) increasing the initial relative Lorentz factor 7q also increases the efficiency. 
However, it also increases the peak energy of the burst. For instance, increasing 7q from 2 to 5 in simulation 
No. 1 increases the total emission by a factor of > 30 i.e. Eiso ~ 10^^ erg, but the peak of energy is at 
> 10 MeV which is in contradiction with both Swift-BAT and Fermi-GBM and Fermi-LAT observations. 
They find the peak energies in the range of few hundred keV for long bursts and ~ 1 — 2 MeV in short hard 
bursts. The shape of the light curve also becomes strange. It includes a very short hard ~ 0.1 sec peak at 
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Figure 18. Evolution of kinetic energy of colliding shells in the rest frame of a far observer at the redshift as source. 
Left: simulation No. 7 with B = 50 kGauss and densities 3 times larger than what is reported in Table 1. This burst 
is very bright with Er 10 x 56 erg i.e. brighter than GRB 080319B and GRB 080607. Right: simulation No. 10, a 
bright short burst. Total kinetic energy (full line), kinetic energy of fast shell (dash-dot), kinetic energy of slow shell 
particle passed through shock front (dash line), magnetic energy inside fast shell assuming uniform density (dash-dot), 
and accumulative emitted energy (dot). The x symbol indicates the value of total emitted synchrotron emission (its 
position on x-axis is arbitrary). Fast variations in magnetic field energy and accumulative energy is due to abrupt 
change in parameters as can be found in Table 1. 



the beginning, then becomes softer and longer peaks later. This setup is not a usual characteristic of GRBs. 
Although in few apparently long bursts such as GRB 060614 [36, 61] very short spikes have been observed, 
they are rare and not a common feature of GRBs. Therefore, the range of parameters that can increase the 
efficiency is quite restricted by other constraints. 

The total kinetic energy of the shock for an observer on the central - for whom shells have spherical symmetry 
is the sum of initial kinetic energy of the fast shell and the kinetic energy of material from slow shell that 
are absorbed by the shock: 
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where Egsc is the total energy of emitted synchrotron self-Compton and T is shock duration until total 
coalition of shells. We use total duration of each simulation as T. But this is somehow arbitrary because it 
does not correspond to total coalition time. Nonetheless, the contribution of slow shell in the total energy 
is small, see Fig. 18 for examples, thus its value does not significantly affect the absolute efficiency, see the 
caption of this figure for more details about these examples. Note that an observer in the center and a far 
observer do not find the same angular distribution for the emission. Due to the relativistic modification 
of emission distribution and different Doppler shift of photons emitted at different angles, the emitted and 
received photons do not have the same angular fiux distribution, see e.g. [103, 113] and one has to take into 
account this difference in the determination of absolute efficiency (4.4). Equation (3.32) for the synchrotron 
power corresponds to the received power, see equations 50 and 56 in [143]. For determining Eggc we cannot 
use equation (3.32), but must use the angular power distribution, transfer it to emitted distribution using 
the relation between emitted and received powers [103]: 

(4.5) 
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where 9' is the angle between photon direction and hne of sight of the observer in the slow shell frame. After 
applying this to equation 56 in [143] and integration over theta' one finds that at first order of approximation 
and when terms proportional to G{r) in (3.32) are negligible: 

(4.6) 

dw 

This relation can be understood if we remind that for far observers emission is beamed in their direction. 
Therefore, the integration for determining total emitted energy must be performed on the beam shape profile 
rather than on a sphere. 

The efficiency of some simulations are reported in Table 2. The absolute efficiency is the ratio of total 
emitted energy calculated above to total kinetic energy from equation (4.3). The table shows also the 
weighted efficiency by the maximum values of eg for each simulation as well as the ratio of emitted energy 
to the kinetic energy of material fallen to the shock front. For presentation in this table, we have chosen 
one example of long bursts and one example of short bursts without and with various values of magnetic 
field. Two other example are presented which show how the variation of other improtant parameters such 
as 7' and eg affects the total emission and its efficiency. Fig. 19 show the spectrum, light curves and energy 
of various components for variety of choices of 7', and duration of regimes. As expected this efficiency is 
usually small, only few percent. Large value of weighted efficiencies, specially when a large external magnetic 
field is present means that the radiative efficiency of shocks is high and most of the energy transferred to 
electrons is emitted very rapidly. As for the absolute efficiency, we should remind that the efficiency of 
synchrotron emission in a shock is always small and this can be understood as the following: Considering 
a neutral plasma, for every accelerated electron that falls to the shock front, one or more baryons with a 
rest mass ~ 2000 times larger than electrons also falls - see Fig. 20 for an examlpe of /?' evolution. Fallen 
baryons make the fast shell heavy, decreases relative Lorentz factor, and thereby fields and synchrotron 
emission. Only the small fraction of kinetic energy which is transferred to fields can be emitted. On the 
other hand, this fraction cannot be high because larger the fields, larger and faster emission, thus faster 
energy dissipation and mixing of two shells which reduces fields and emission. Therefore, this process is 
self-regulatory. Despite low efficiency of emission the energy budget of simulations in Table 2 show that 
even for brightest bursts, the total kinetic energy of the outflow is of the order of energy released in normal 
supernovae, thus plausible. 

Apriori the end of the shock is when the two shells are completely mixed and move together with the same 
Lorentz factor. In reality this complete mixing occurs much later than the end of high energy emission. This 
can be seen in Fig. 18. It shows that during the evolution of the burst most of the energy is emitted at early 
times - consistent with observed FRED-like peaks. We also note that in absence of an external magnetic 
field emission efficiency is usually very low. These features can be described as the following: Initially, the 
discontinuity of density and other physical characteristics of the shells is large and turbulence important. 
Formation of the shocked and active regions, and thereby energy dissipation changes the discontinuity to a 
gradual variation that connects high and low density regions. At this point, the lack of sufficient turbulence 
stops high energy emission and only low energy synchrotron emission due to a weak remnant magnetic or 
due to the weakening external field would persist. Evidently our simple phenomenological model does not 
include these details. Nonetheless, parameters of this model are meant to present these processes, thus 
we expect that their evolution grossly reproduce missing details. Fig. 20 shows an example of evolution 
of various quantities in the shocked region. The peak of emission is during the time interval around the 
maximum of UJm- 

If the width of the fast shell is much shorter than slow one, it dissipates completely, but this can take a very 
long time. If the fiow of fast material continues, accumulation and compression inside and in front of the 
transition region form another shock front and a new spike of high energy emission. Therefore, an emission 
gap is expected from the formation of a gradual transition region to accumulation of material and formation 
of a new shock. During this period the fast shell propagates freely and without making any radiation. These 
qualitative descriptions along with our simulations suggest a solution for the efficiency of energy transfer 
to electrons and thereby the synchrotron emission. Because shocks are efficient only for a short time, if 
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Figure 19. Simulation No. 6 with various values for 7' and r/ro: 7' = 2 (top), 7' = 4 (middle). 7' = 4 and for 4 
regimes define in Table 1 r/rg = 1.5 — 5 — 5 — 5. Other parameters are the same as what is reported in Table 1. Note 
that light curves are presented for an observer at the redshift of the burst. Highest efhciency is obtained with 7' = 2, 
and highest energy output for 7' = 4 and r/rg ~ 1.5 — 5 — 5 — 5. Note that variation oi r/ro without changing indices 
of €e and €b varies the peak value of the latter and thereby the synchrotron emission. We have also tried 7' = 20 and 
7/ = 50. This leads to the same initial F, but synchrotron emission has a low flux, is soft except for a very short peak 
at the beginning, and efficiency is low ^ 0.1%. For these reasons they are not shown here. 



anisotropies are strong at small spatial scales but are very frequent, there is a large probability that the 
fast shell collides with another shell before significant suppression of discontinuity. In this way the high 
efficiency of the emission can persist for much longer time and a larger fraction of the kinetic energy can be 
transferred to emission. This means that long bursts consist of overlaping short bursts. 



5 Outline 



In this work we introduced an external magnetic field in a previously developed formulation for relativistic 
shocks and synchrotron self-Compton emission as a model for the prompt emission of GRBs. Simulation of 
a number of bursts with properties similar to observed ones showed that due to the fast variation of many 
quantities in this model, oscillations of an external field do not leave significant evidence of their presence in 
the PDS of light curves. Therefore, the lack of a clear oscillatory signature is not the evidence for absence of 
an external magnetic in GRBs. In fact according to the present knowledge of the physics of massive stars - 
presumed progenitors of GRBs - some sort of magnetism must exists in their environment. We also argued 
that a distance of order 10^*^"^^ cm for the shock from the central engine produces more realistic bursts and 
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Figure 20. Evolution of quantities determining dynamics for simulation No. 10. All quantities except F and Ar are 
in slow shell frame. The time is in the rest frame of an observer at redshift of the burst. 



solves the problem of the undetected photospheric emission of a standard fireball. 

The spectra of simulated bursts show that the spectral properties of observed GRBs can be explained by 
internal shock synchrotron self-Compton model. In particular, diversities of the high energy wing of the 
spectrum and flatness of its low energy wing seem to be related to the properties of electrons Lorentz 
factor distribution. The spectrum of very hard and usually short bursts can be explained with a power-law 
distribution which continues to very large energies without any cutoff. Softer and usually longer bursts 
have either broken power-law or power-law with exponential cutoff. The shallow slope of observed GRBs 
is also most probably related to the distribution of electrons and its evolution. It can be explained with a 
distribution that continues smoothly to very low accelerations rather than having an abrupt break as usually 
considered. 

Inverse Compton scattering can have an observable contribution at low energies, consistent with the results 
of [22], but at high energies its contribution is orders of magnitude smaller than synchrotron. Nonetheless, 
it is not straight forward to estimate the total contribution of inverse Compton scattering to GRBs because 
there can be more low energy photons of non-shock origin in their environment. They can increase the flux 
of inverse Compton scattered photons. 

As for the reason for much longer emission in high energies we suggested that it should be related to a 
resonant movement of highly accelerated electrons with the electromagnetic standing wave created by the 
shock. This coherent movement delays the time of passage from an acceleration zone to a dissipation zone, 
increases the fly time of high energy electrons and postpones their dissipation. 
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A Analytical approximation of 7V( 
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When r/i/(3 — k) is large we can formally expand the integrand of the integral in (3.24) and integrate term 
by term: 



AA(o)(r') 



Ar' 



{3-k) 



/3, 



/T-2+- 



/3; 



(0) 



11 



2FI 



iFii 





T-2 + 




T-1 + 




3 — K, 


3 — K 


3 — K 


Vi 


T-2 + 


m 


T-1 + 


Vi 


3 — K 


3 — K 


3 — K 



/3, 



° )] 



(3 - k)V 



+ 



(i-i) 



/3(o) 



iFi[ 



3 — K 



3 — K 



^ + 75 ■,T + l + - ;1 

3 — K 3 — K 



/3o 



{3-k)V' 



3 — K 



T + 1 + 



Vi ^(0)'- 



3-k' 13' 



(3 - k)V 



)] 



+ ... 



(A.l) 



where /3q in the initial (3 = v/c of the fast shell with rest to the slow shell and /3|q)(7'') is the zero order 
solution i.e. when the energy loss due to synchrotron emission is ignored: 
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2F1 is the hypergeometric function and can be expanded as a polynomial of its last argument. When the 
latter is less than 1 and the power of the terms in the polynomial are positive, they converge rapidly to zero. 
When the last argument in 2-^1 is larger than one an analytical extension of this function with negative 
power in the polynomial expansion exists. Therefore, the dominant term is always /3^q^ in front of each 2F1 

term. For small rji this approximation is not valid. In this case we expand (1 + ( (3_k)p ~ 
obtain: 
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B Projection of precessing field on the shock front surface 

Figure 21-a presents the geometry of the magnetic field with respect to the shock front, ip is the angle 
between the precession axis of the external magnetic field and the shock front surface S. The opening angle 
of precession is p. For the sake of simplicity we first determine the projection of precession on a surface 



-44- 




Figure 21. a) Precessing external magnetic field and its projection on S the shock front source - surface orthogonal 
to the propagation direction. The magnetic field B precesses with period w around a constant axis having angle ^ 
with E. p is the precession angle and B± is the projection of B on Sigma, b) The precession loci of B±^ on Sigma. 



at rest with respect to the central engine and parallel to the shock front . Then, we transform projected 
quantities to the slow shell frame. 

Projection of the circular precession orbit of the field is an ellipse with half diagonal a and b where a is 
parallel to the intersection between shock front and rotation surface: 

a = Bsinp b = B sin ip sin p (B.l) 

Figure 21-b shows the projection of the magnetic field and its component parallel to the surface of rotation 
on S. Using the ellipse equation we find following expressions for x and y components of the field parallel 
to rotation surface and the distance between points A and C (see Fig. 21): 

X = \B\sin pcos{uj{t — to)) y = \B\sin psinip sin{uj{t — to)) AC = \B\ cos{^ + p) (B.2) 

Finally: 

Bl = {AC^ + b + y)^ + x^ = \B 



cos(V' + p) + sinpsinV'l 1 + sin(tj(t — to)) 



(B.3) 



It is assumed that no large scale electric field is generated in the engine, thus the transformation of magnetic 
field to a frame moving in the direction of z-axis with a Lorentz factor of F is very simple [47]: 

B', = B, B'^ = TB^ (B.4) 

In contrast to kinematic quantities, the magnetic field along the boost direction is preserved, but the 
transversed field is boosted. Schematically, the observer in the boosted frame see that charge particles in 
the engine are moving with a Lorentz factor of F creating a current which in its turn creates a magnetic 
field in the surface transverse to current direction. 



C Plasma instabilities in presence of a background magnetic field 



In the simple simulations presented in Sec. 4 it was only possible to consider the effect of an external 
magnetic field on the synchrotron emission. However, its influence on the development of instabilities in the 

^''For simplicity here we identify the active region with the shock front because the latter is effectively very close to a flat 
plane and assume that the active region is thin. The case of a thick active region can be considered as superposition of multiple 
thin emitting layers. 



plasma, the properties of instabilities and their evolution are much more complex. Unfortunately, it is very 
difficult to formulate these processes and find analytical or even numerical solutions for various physical 
quantities. For this reason in this section we briefiy and qualitatively review the properties of a plasma in 
presence of a magnetic field according to latest PIC simulations. As mentioned before, due to computational 
restrictions, none of these simulations fulfills all the requirements necessary for the production of a GRB. 
But each of them targets some aspects of the physics of ultra-relativistic shocks. Conclusions of these studies 
are important for understanding the reasons for absent or very weak signature of a magnetic field in GRB 
data and are complementary to the conclusions obtained from simulations presented in this work. 

In absence of an external magnetic field in a plasma environment, Weibel instabilities [130] produce a co- 
herent field only at scales much shorter than the size of a shell. In the GRB context, the shock between 
relativistic shells can be considered as a beam-background plasma process [5, 64]. Anisotropics in the distri- 
bution of oppositely charged particles across the shock front induces a magnetic field that grows for some of 
fluctuation modes k until particles get trapped in this field and their growth stops. The growth rate of 
the Weibel instabilities is defined as the imaginary part of the wave number = '^muj{k). The amplitude 
of the induced magnetic field depends on which in its turn depends on the distribution of both ions 
and electrons in the two sides of the shock front [132]. For thermal distributions an explicit solution for F^ 
exists [112, 134]. For relativistic shocks in absence of a large scale magnetic field F^^'^^'* ~ ^'pe/l'^^'^ where 7' 
is the relative Lorentz factor of the shells and oj^g = (47re^?ig/me)^/^ is the plasma frequency of electrons [64]. 
For typical densities used in our simulations, the time scale for the growth of Weibel instabilities is very 
short ~ 10~^ sec for a far observer. 

In presence of an external magnetic field the dispersion relation of the plasma is more complex and depends 
on the angle between the external field and fiuctuation modes. In fact due to the magnetic force particles 
of the slow shell penetrate only a short distance into the fast shell region before they are defiected [51]. 
This makes the mixing of the two shells more difficult, but increases the efficiency of synchrotron emission. 
Moreover, the external magnetic field induces a phase shift between electric and magnetic potentials that 
helps to separate the acceleration electrons from their dissipation. These results are also confirmed by 
3D-simulations [72]. 

Startsev, et al. [112] have studied a special case in which the external field is parallel to the direction 
of the movement of the plasma, both analytically and numerically. They show that depending on 9 the 
angle between plasma wavenumber k and the external magnetic field as well as the relative importance of 
parallel and transverse modes, instabilities can be either Weibel-like, electrostatics lower-hybrid or higher- 
hybrid. Detailed definitions and properties of these instabilities can be found in [112]. Weibel instabilities 
occur when 6 is very small and only for long parallel modes (parallel to the field and to the bulk velocity 
of plasma) and A^y <^ k±. For larger k^^ and the same 9 other type of instabilities are formed. The main 
physical difference between these instabilities is whether ions (protons) of fast shell or slow shell stream with 
respect to electrons and create the induced fields. Notably, although in presence of an external magnetic 
field instabilities continue to grow, their rate saturates - see Figure 4 in [112]. The saturated growth rate 
(width) for P' J^J uj'^^ < 1: 

F_ = -^^^^ (C.l) 

Here uj'^^ = eBext/nie is the cyclotron (Larmor) frequency of electrons in the slow shell, uipf. = (47re^ne/rfie)^''^ 
and ^'pef are plasma frequency of protons in the slow and fast shells respectively. We assume a plasma 
composed of protons and electrons. For the range of densities and magnetic fields we have used in the 
simulations, the plasma of slow and fast shells are in the saturated regime. This means that the shock 
produce a transverse magnetic field in the active region, and Fermi or other similar processes [115] accelerate 
electrons and produce synchrotron radiation, but this field is limited. For simulations in this work the growth 
width F^oo ^ 10~^ eV in slow shell rest frame, corresponding to a time scale of < 10~^ sec in the rest frame 
of a far observer. It is longer than unmagnetized plasma but very short with respect to time resolution of 
available observations. We should also mention that in a magnetized plasma Weibel instabilities are not 
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always dominant [112]. 



Propagation of shock waves and instabilities and maximum achievable induced field depend also on the value 
of ojce/^pe and on Vgh/vA where Vsh is the velocity of the shock front /fast shell and va = B / (ATrnnip)^/'^ 
is the Alfven velocity of the slow shell [51]. For mildly relativistic internal shocks and values of |-B| and 
density of the plasma considered in our simulations, Vsh/vA ^ 1 thus the propagation distance of ions inside 
the shocked region is large. In fact, both simulations and observations find jets and shocks in sources with 
much larger magnetic fields than what is used in these simulations, such as in neutron stars and pulsars, see 
e.g. [49? ]. 

The fact that instabilities become active at distinct range of plasma fluctuations and external magnetic fields 
strength is very important for GRBs. As we mentioned above plasma modes have a quasi-particle behaviour 
with a non-linear and sophisticated relation between frequency and wavenumber. During propagation of the 
ejecta and collision of shells some modes can get a lag/lead with respect to others. Moreover, the amplitude 
of the induced magnetic field depends in a non-linear manner on the external magnetic field which can be 
time dependent. Here we should remind that the wavelength of plasma waves expected to be of order of 
few centimeters for electrons and few thousands of centimeters for protons/ions. Therefore, they can apriori 
explain very short time scale variations - the jitter - of GRBs light curves [63], but cannot explain a coherent 
oscillation or even random variations of multiple-second scale in long bursts. Therefore, these features must 
be due to large scale variations in the properties of the ejecta, e.g. density, Lorentz factor, etc. and/or 
multiple separate ejection of material by the central engines. If the external magnetic field in the direction 
of the ejecta oscillates - for instance due to the precession of the rotation axis of the central engine - the 
induced magnetic field and thereby the synchrotron emission of accelerated electrons vary cyclically. This 
also induces an oscillation into distributions of electrons and ions because dispersion relations depend on 
the external magnetic field. 

Due to the nonlinear relation between an external magnetic field and the induced magnetic field by insta- 
bilities, their oscillation frequencies are not necessarily the same. Moreover, during the propagation and 
collision of shells the nature of dominant instabilities and thereby their dependence on the external magnetic 
field can change. Therefore, one expects that effects such as: dominance of higher harmonics, lags due to 
the differential dispersion relation of modes, intrinsic variation of plasma density along the line of sight, 
etc. smear the coherent behaviour in a large extent. Thus, the probability of observation of a significant 
oscillation in the light curve of GRBs is very small. Due to these complexities, a consistent treatment of 
an oscillating external magnetic field on a GRB emission needs a full plasma physics treatment of these 
processes as well as a detailed knowledge of the jet/ejecta properties which are not available. Thus, for 
studying the effect of an oscillating external magnetic field on the synchrotron emission we have simply 
assumed that the total magnetic field consists of a magnetic field which depends on the properties of the 
plasma and an external oscillating component which can also varies with the distance of the ejecta from 
the central engine. In addition, for further simplification of our calculation we consider the oscillation as 
semi-stationary i.e. we ignore terms that depend on the derivative of the oscillations. 



D Electron cooling 

The usual definition of "slow cooling" and "fact cooling" suggested first in [106] and have been used 
intensively in GRB literature cannot be applied to this model because quantities such as Lorentz factor F 
and magnetic B, to which the Lorentz factor of cooled electrons 7c depends, vary with time. In this case 
the definition of 7c can be changed, for instance to: 

7c = / dtP(7c)/Fmec2 (D.l) 
Jo 

Although the frequency dependence of 7c defined here is apriori the same as what is discussed in [106] , due 
to time-dependence of power P, its variation with time is not the same as what is obtained in [106] . Slow 
and fast cooling regimes are defined by comparing 7c wjth 7^ which is also time-dependent, see Fig. 20. The 



presence of an external magnetic field and electron distributions other than power-law make the situation 
even more different from the simple model used to predict broad band spectrum of GRBs. Moreover, both 
the definition (D.l) and the original definition of 7c depend on unknown and somehow arbitrarily defined 
dynamical time of emission. Therefore, 7c, its comparison with 7^, and classification of processes as fast and 
slow cooling are useful for getting an insight into behaviour of the spectrum expected from a synchrotron 
model but not for analyzing a real burst or a more complex model such as the one presented in this work. 
For these reasons we have not used this classification in our analysis. 

Qualitatively one expects that at the beginning when fields are not yet very strong i.e. eg and are small, 
electrons are in slow cooling according to the convention of [96]. When the amplitude of magnetic field, 
acceleration of electrons, and thereby 7^ increase, electrons distribution becomes harder and they move to 
a fast cooling regime. Finally when the shock becomes weak, distribution of electrons becomes softer, they 
go back to a slow cooling regime, and emit dominantly at low energies. See Fig. 20 where an example of 
the evolution of quantities determining characteristics of the shock and synchrotron emission are shown. 

In Figs. 13 and 14 spectra of different time intervals for some simulated bursts are plotted separately. They 
show the gradual movement of the peak to lower energies. This is similar to switching from "fast cooling" 
to "slow cooling". In fact the only difference between fast and slow cooling regimes - for simple power-law 
distribution - is the slope of the flux around the peak. For slowly varying afterglows which was the subject 
of study in [96], fast and slow cooling regimes occur in well separated time intervals. In the case of fast 
varying prompt emission both regimes occur during the short duration of prompt emission because the 
distribution of electrons changes with time (similar to our simulations). Therefore, the spectrum is formed 
from superposition of multiple fast and slow cooling regimes, and population of accelerated and ones which 
have kept more or less their original distribution - thermal or non-thermal. It is the reason for the absence 
of features corresponding clearly to one or the other regime in observed spectra. 
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Table 1. Parameter set of the simulated models. 
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lOphase = wio in (B.2) 

iiSome plots shown in the text differs only in one variable. In this case the model 
I2parameters is shown in this table only once. 

l3Each line corresponds to one simulated regime. Horizontal lines separate independent simulations, dashes mean either the quantity 
is evolved from its initial value according to evolution equations ^gpcribed in the text or kept constant because its variation induces 
discontinuity to light curves and the model is not adapted for their variation. Parameters t and 5 arc used in different Ar' evolution 
models and are mutually exclusive. 
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^^This simulation has a 7' larger than others and T f much smaller. Apriory it must have higher efficiency and indeed 
it has, but its efficiency yet only few percent and only 56% largeer than simulation No. 7 with same magnetic field. 
^^The idea behind this simulation is to see whether large eg and 7' but flat electron distribution which reduces 7^, 
see (3.10) and (3.12) produces a higher efficiency. It seems that this is not the case. The peak energy spectral shape 
are also inconsistent with observations. 
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